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Abstract 


'Alin  a  combination  tutorial  and  research  paper,  techniques  for 
combating  the  effects  of  multipath  on  high-rate  data  transmissions 
via  radio  are  explored.  The  tutorial  aspect  of  the  paper  presents: 
(1)  a  heuristic  outline  of  twenty-five  years  of  work  on  the  theory  of 
anti-multipath  radio  receivers  and  (2)  a  summary  of  a  statistical 
model  of  urban/suburban  multipath.  The  research  section  of  the  paper 
presents  results  of  analyses  and  simulations  of  various  candidate 
receivers  indicated  by  the  theory,  as  they  perform  through  urban/ 


suburban  multipath.  A  major  result  shows  that  megabps  rates  through 


urban  multipath  (which  typically  lasts  up  to  5  ps)are  quite  feasible 


r 
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I.  Introduction 

Some  digital  radio  systems  must  operate  through  an  extremely  harsh 

multipath  environment,  in  which  the  duration  of  the  multipath  exceeds 
* 

the  baud  length.  Two  disciplines  combine  to  shed  light  on  receiver 
design  for  this  environment:  the  modeling  and  simulation  of  multipath 
channels  and  the  theory  of  multipath  and  other  diversity  receivers. 

/  In  this  paper,  we  first  present  a  tutorial  review  of  pertinent 
aspects  of  both  underlying  disciplines.  We  then  carry  out  rough  but 
indicative  analyses  of  the  performances  of  two  promising  binary  anti¬ 
multipath  systems.  Finally  —  since  the  analyses  contain  a  number  of 
oversimplifications  that  make  them  heuristic  rather  than  definitive  — 
we  present  results  of  computer  simulations  of  the  two  proposed  con¬ 
figurations  and  others,  as  they  operate  through  simulated  urban/suburban 
multipath. 

A  heavy  debt  of  gratitude  is  due  the  National  Science  Foundation 
for  its  support  of  underlying  basic  research  on  multipath  channel 
characterization  and  simulation  at  the  University  of  California  (Berkeley). 
Analyses  and  simulations  of  particular  systems  were  supported  by  the 
Advanced  Research  Projects  Agency  (ARPA)  through  SRI  International. 

Finally,  the  author  acknowledges  with  pleasure  his  students  U.  Charash, 

H.  Hashemi,  M.  Kamil  and  H.  Suzuki,  upon  whose  Ph.D.  research  this 
paper  heavily  relies. 


An  example  is  the  ARPA  Packet  Radio  network  currently  being  designed 
and  constructed. 
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1 1 .  Modeling  Multipath  l’rnpng.u  1qi> 

Ultimately,  a  reliable  multipath  model  must  be  based  on  empirical 
data  rather  than  on  mathematical  axioms.  Two  types  of  data  are  available. 
The  more  common  type  give  the  results  of  narrowband  or  CW  measurements, 
in  which  only  a  single  fluctuating  variable  —  a  resultant  signal  strength 
—  is  measured  [11].  Although  the  fluctuation  of  this  strength  variable 
depends  on  reception  via  multiple  paths,  these  paths  are  not  resolved  by 
the  measurements.  I'**  shall  denote  the  results  of  such  measurements  a9 
"fading"  data  rather  than  multipath  data,  because  they  determine  a 
fading  distribution  of  the  single  strength  variable  e.g.,  Rayleigh 
lognormal,  Rice,  etc. 

Wideband  experimental  data  that  characterize  individual  paths  are 
less  common  [5],  [10],  [ 191,  [33],  [37].  In  order  to  resolve  two  paths 
in  such  measurements,  the  sounding  signal's  bandwidth  must  he  larger  than 
the  reciprocal  of  the  difference  between  the  paths '  delay.*:.  Although 
bandwidths  of  100  MHz  or  more  have  been  used  in  exceptional  circumstances 
to  resolve  path  delay  differences  of  less  than  10  ns  [10],  the  bulk  of 
available  data  derive  from  10-MHz  bandwidths  or  less  [5],  [19],  [331, 

[37].  In  the  latter  measurements,  paths  separated  by  delays  of  more  than 
100  ns  are  resolved;  multiple  paths  with  smaller  separations  are  seen 
as  single  paths. 

The  nature  of  the  multipath  measurements  depends  somewhat  on  the  use 
envisoned  for  them.  If  understanding  of  the  effect  of  the  multipath 
channel  on  CW  transmissions  is  required,  measurements  that  show  Doppler 
effects  may  be  important  [5],  and  those  are  reasonably  related  to  a 
scattering-medium  model  of  the  channel  [1],  [12].  For  high-rate  packetized- 
data  transmission,  for  vehicle-location  sensing,  and  for  other  "bursty" 
transmissions,  measurement  of  sequences  of  "impulse  responses"  of  the 
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propagation  medium  suffices. 

The  simulations  of  data  reception  that  are  presented  in  a  later 
section  are  based  on  the  "impulse  response"  approach,  and  it  is  to  this 
type  of  model  that  we  restrict  ourselves.  In  order  that  the  model  and 
the  simulations  themselves  be  fully  understood,  we  shall  review  here  the 
experiments  underlying  the  model.  These  were  performed  in  urban/suburban 
areas  [32],  [33]. 

A .  The  Underlying  Experiment s 

Pulse  transmitters  were  placed  at  fixed,  elevated  sites  in  the 
San  Francisco  Bay  Area.  Once  per  second,  these  would  simultaneously 
send  out  100-ns  pulses  of  carrier  at  488,  1280  and  2920  MHz.  The  pulses 
were  received  in  a  mobile  van  that  moved  through  typical  urban/suburhan 
areas,  recording  on  a  multitrace  oscilloscope  the  logarithmically  scaled 
output  of  the  receiver's  envelope  detectors  (see  Fig.  1).  Since  the 
oscilloscope  was  triggered  by  a  Rubidium  frequency  standard  that  was 
synchronized  with  a  similar  unit  at  the  transmitters  prior  to  each 
experimental  run,  absolute  propagation  delays  could  be  measured  within 
experimental  accuracies  of  better  than  20  ns. 

Four  series  of  experiments  were  performed  in  the  following  typical 
urban/suburban  areas: 

A:  Dense  high-rise  —  San  Francisco  financial  district 

B:  Sparse  high-rise  —  Downtown  Oakland 

C:  Low-rise  --  Downtown  Berkeley 

D:  Suburban  —  Residential  Berkele.y 

In  each  area,  regions  of  dimensions  roughly  300-1000  feet  (along  the 
transmitter-receiver  line;  of  sight)  by  2300-4000  feet  (tangential  to  line 
of  sight)  were  exhaustively  canvassed,  with  care  taken  to  include  proper 
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LOS  DELAY  8EYOND  LINE  OF  SIGHT  (LOS)  IN  MICROSECONDS 


Figure  1.  Example  of  measured  multipath  profiles  for  a  dense  high 

rise  topography.  (a)  Top  to  bottom:  2920,  1280,  488  MHz. 
Vertical  scale:  35  dB/cm.  Horizontal  scale:  1  ps/cm. 
Different  apparent  LOS  delays  are  due  to  difference  in 
equipment  delays.  (b)  Middle  trace  of  (a)  on  a  linear 
scale. 


topographic  cross-sections:  intersections,  mid-blocks,  points  at  which 
the  transmitter  site  was  visible  or  occluded,  etc.  About  one  thousand 
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frames  of  data  of  the  type  shown  in  Fig.  1(a)  were  obtained  in  each  area. 
B.  A  Fundamental  Model 

The  model  upon  which  data  reduction  was  based  was  one  first  posed 
in  [27].  In  this  model,  it  is  assumed  that  a  transmission  of  the  form 

jw  t 

s(t)  =  Re [ o ( t ) e  °  ]  (1) 

will  be  received  as 

jw  t 

r (t)  =  Re[p (t)e  °  ]  +  n(t)  (2) 

where 

P(t)  =  V  a  o(t-t.)  e  (3) 

In  (l)-(3),  o(t)  is  the  complex  envelope  of  the  transmission,  i.e.,  |o(t) 
is  its  amplitude  modulation  and  tan  ^[Im  o(t)/Re  o(t)]  is  its  phase 
modulation.  The  transmission  is  received  via  a  number  K  of  paths,  the 
kth  of  which  is  characterized  by  three  variables:  its  strength  a^,  its 
modulation  delay  t^»  and  its  carrier  phase  shift  0^.  The  waveform  n(t) 
is  an  additive  noise  component. 

It  is  assumed  that  all  paths  are  resolvable,  i.e.,  they  are  counted 
as  different  paths  only  if  the  condition 

It^-t^j  >  1/W  for  all  k  ^  l  (4) 

holds,  where  W  is  the  transmission  bandwidth.  Distinct  paths  in  the 
physical  medium  that  violate  this  resolvability  condition  are  not  counted 
separately,  since  they  cannot  be  distinguished  by  a  measurement  of 
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bandwidth  W.  Instead,  any  two  paths  —  call  them  k^  and  —  for  which 


jt.  -t  I  <  1/W  are  considered  as  a  single  path  in  (3),  with  a  common 
1  2 

delay  t,  a  t,  5  t,  and  a  strength/phase  combination  given  by 

K  lC  -  Nfl 


JO 


'1 

j0, 
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k  A 


\  e 
K1 


+  a.  e 
k2 


It  is  the  triplet  {t^a^t0^  that  is  to  be  determined  for  each 
"resolvable"  path.  To  be  sure,  if  a  continuum  of  paths  existed,  it  would 
be  difficult  uniquely  to  cluster  the  "subpaths"  into  paths.  But  many 
media,  including  the  urban/suburban  one,  have  a  natural  clustering 
—  e.g.,  groups  of  facades  on  buildings  —  that  make  the  model  feasible. 


C.  A  Discrete-Time  Approximation  to  the  Model 

In  addition  to  the  additive  random  noise  n(t)  in  (2),  the  received 

K-l 

signal  r(t)  is  therefore  characterized  by  the  random  variables  (t^Q  » 

K-l  K-l 

^ak^0  ’  ^k^O  and  K‘  Tke  Purpose  data  reduction  from  the  "multipath 

profiles"  exemplified  by  Fig.  1  was  to  obtain  statistics  of  these  random 
variables  upon  which  to  base  a  simulation  program.  A  number  of 
generations  of  statistical  models  —  based  both  on  the  data  and  on 
physical  reasoning  when  the  data  were  insufficient  or  undecisive  —  ensued 
[32],  [25],  [26],  [8],  [9].  The  following  final  version  emerged. 

Each  multipath  profile  starts  with  the  line-of-sight  (LOS)  delay,  which 
is  chosen  as  the  delay  origin.  Since  the  resolution  of  the  original 
experiment  is  100  ns,  the  delay  axis  is  made  discrete  by  dividing  it  into 
100-ns  bins,  numbered  from  0  to  70.  Bin  0  is  centered  on  LOS  delay, 
subsequent  bins  being  centered  on  multiples  of  100  ns.  The  delay  of  any 
physical  path  lying  in  bin  f  is  quantized  to  100C  ns,  the  delay  of  the  bin's 
center.  Figure  2  shows  the  bin  structure,  a  multipath  profile,  and  the 
resulting  discrete-time  path-delay  structure.  Notice  that  only  paths 
with  delays  less  than  7.05  us  beyond  LOS  delay  are  encompassed  in  this 
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Figure  2.  Discrete-time  model.  (a)  Division  of  excess  delay 
axis  into  71  100-ns  bins.  (b)  A  typical  multipath 
profile.  (c)  Discrete-time  path-delay  indicator 

string, 
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model;  experimental  evidence  shows  that  significant  paths  with  larger 
delays  are  highly  improbable. 

K—  i 

The  path-delay  sequence  (t^^  is  approximated  by  a  string 

{t^l^of  0's  and  l's,  as  shown  in  Figure  2(c).  If  a  path  exists  in  bin 

l,  t  =  1;  otherwise  b  0.  In  the  sample  string  in  Figure  2(c),  only 

Trt,T_,T,  ,t  . . t,„  ,x  are  nonzero,  corresponding  to  the  quantized  path 

U  J  4  o  oz  b/ 

delays  t^  =  0,  t^  =  300  ns,  "  400  ns,  t^  =  600  ns,...,t^_2  =  6300  ns, 
tR  ^  =  6700  ns.  (t^  is  the  value  of  t^,  as  quantized  to  the  nearest 
100  ns. ) 

Associated  with  each  nonzero  is  the  corresponding  pair. 

Thus,  the  discrete-time  model  is  completed  by  appending  to  the  string 

K-l 

a  set  of  strength-phase  pairs  {  (a^,  6^)  where  the  index  k  refers  to 

the  kth  nonzero  entry  in  the  string.  This  is  shown  in  Figure  2(b). 

The  discrete-time  model  of  Figure  2  pictures  the  multipath  profile 
at  a  single  point  in  space.  A  sequence  of  such  profiles  is  needed  to 
depict  the  progression  of  multipath  responses  that  would  be  encountered 
by  a  vehicle  following  a  track  such  as  shown  in  Figure  3(a).  One 
imagines  points  1, 2 , . . . ,n, . . . ,  arbitrarily  placed  on  the  track,  at  each 
of  whicl  a  multipath  response  is  seen.  The  discrete-time  versions 
of  these  responses  are  arrayed  in  Figure  3(b),  where  an  additional  spatial 
index  n  has  been  superscribed  on  all  variables. 

One  begins  to  recognize  the  complexity  of  the  model  and  of  the  required 
reduction  of  experimental  data  on  realizing  that  in  addition  to  the 
need  for  first-order  statistics  of  the  random  variables  a^n^ ,  0^°^. 

and  (where  0  <_  l  70;  0  <_  k  <_  K^-l;  1  <_  n  <  ®) ,  there  are  two 

dimensions  along  which  at  least  second-order  statistics  are  necessary: 
temporal  and  spatial.  For  each  profile  (fixed  n) ,  there  are  temporal 
correlations  of  the  delays,  strengths  and  phases  of  the  several  paths; 
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Figure  3.  Spatial  variation  of  the  multipath  profile.  (a)  A 

vehicle’s  track,  with  spatial  sample  points.  (b)  The 
sequence  of  multipath  profiles  at  the  sample  points 
on  the  track. 
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In  addition, there  are  spatial  correlations  of  these  variables  at 
neighboring  geographical  points. 

The  reduction  of  experimental  data  f  33  ] ,  [25],  |9|  and  physical 
reasoning  led  to  the  following  model,  wh ii  h  was  the  basis  for  simulation. 

(1)  The  string  of  the  nth  profile  is  a  modified  Bernoulli 

sequence,  in  which  the  probability  of  a  1  in  the  it h  place  depends 

on:  (a)  the  value  of  l ;  (b)  whether  a  1  or  a  0  occurred  in  the 

s  t 

(£-1) —  place  of  the  same  profile;  (c)  whether  a  1  or  a  0  occurred 

s  t 

in  the  ith  place  of  the  (n-1)- —  profile. 

(2)  The  strength  a^^  of  the  kth  path  of  the  nth  profile  is 

* 

conditionally  log-normally  distributed,  the  conditions  being  the 

s  t 

values  of  strengths  of  the  (k-1)-—  path  of  the  nth  profile  and 

s  t 

of  the  path  with  the  closest  delay  in  the  (n-l)- —  profile;  appropriate 
empirically  determined  correlation  coefficients  govern  the  influence 

these  conditions  exert  on  n,^.  The  mean  and  variance  of  the 

k 

distribution  of  a^^  are  also  random  variables,  drawn  from  a  spatial 
random  process  that  reflects  large-scale  inhomogeneities  in  the 
multipath  profile  as  the  vehicle  moves  over  large  areas. 

(3)  The  phase  0^°^  of  the  jcth  path  of  the  nth  profile  is  independent 

of  phases  of  other  paths  in  the  same  profile,  but  has  a  distribution 
depending  on  the  phase  of  a  path  with  the  same  delay  in  profile  (n-1), 
if  there  is  such  a  path;  if  no  such  path  exists,  is  uniformly 

distributed  over  (0,210. 


Actually,  Suzuki  [25]  showed  that  paths  with  small  delays  (beyond  LOS 
delay)  were  better  modeled  by  Nakagami  distributions,  but  Hashemi  [9] 
was  forced  to  approximate  these  log-normally  because  of  the  complexity 
of  the  simulation  program. 
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(4)  The  spatial  correlation  distances  of  the  variables  just 
described  vary  considerably,  ranging  from  less  than  a  wavelength 
for  the  O^'s,  through  tens  of  wavelengths  for  the  a^'s  and  t^'s, 
to  hundreds  of  wavelengths  for  the  means  and  variances  of  the  a^'s. 

These  statistics  are  more  fully  explained  in  [8],  [9]. 

D.  The  Simulation  Program 

Hashemi's  simulation  program  SURP,  based  on  the  statistics  just 
outlined,  generates  sequences  of  multipath  profiles,  as  depicted  in 
Figure  3.  If  one  were  to  examine  a  sequence  of  such  profiles,  he  would 
see  paths  appearing  and  disappearing  at  a  rate  depending  on  the  spacing 
of  points  on  the  vehicle's  track  (Fig.  3(a)).  Profiles  at  only  slightly 
separated  ge’ographical  points  would  look  very  similar,  with  high 
correlations  of  path  delays  and  strengths  (and,  for  very  close  points, 
phases).  Profiles  at  greatly  separated  points  would  not  only  have  grossly 
dissimilar  { t ^ } ,  {a^.S^}  strings,  but  the  gross  strength  statistics  of 
these  strings  (e.g.,  the  average  strength  of  the  paths  in  a  string)  would 
be  dissimilar,  reflecting  the  spatial  inhomogeneity  incorporated  into  the 
model.  The  "motion  picture"  of  simulated  profiles  just  described  is  in 
fact  very  much  like  experimental  data  [37 ] • 

The  simulation  program  can  be  run,  using  empirically  determined 
parameters,  for  each  of  the  three  frequencies  and  four  areas  of  the 
original  experiment.  Long  sequences  of  strings  were  generated  for 
each  of  the  twelve  f requency/area  combinations,  assuming  that  the  points 
on  the  vehicle  track  are  uniformly  spaced  by  distance  d.  An  example  of 
such  sequences  is  given  later  in  Fig.  25.  For  various  values  of  d,  the 
statistics  of  the  simulated  sequences  were  then  compared  with  the  original 
empirical  statistics.  F.xcellent  agreement  was  obtained  [8],  [9].  (See 


Figure  4  for  examples.) 


BUT  1  DNS 


probability  that  the  strength  of  a  path 
in  the  indicated  delay  interval  is  less 
than  the  abscissa  value. 


Sequences  generated  by  SURP  provide  a  data  base  with  which  to 
perform  experiments  with  urban/suburban  radio  systems.  It  should  be 
noted  that  initial  simulation  experiments  on  urban/suburban  radio 
ranging  and  location  systems,  using  a  very  rudimentary  propagation 
simulation  program,  gave  results  which  compared  extraordinarily  accurately 
with  actual  hardware  experiments  [34],  This  initial  success  encouraged 
the  development  of  the  more  elaborate  simulation  capability  just  described 


III.  Design  of  Multipath  Receivers 

Multipath  reception  is  one  form  of  diversity  reception,  in  which 
information  flows  from  transmitter  to  receiver  via  the  natural  diversity 
of  multiple  paths  rather  than  via  the  planned  diversity  of  multiple 
frequency  channels,  multiple  antennas,  multiple  time  slots,  etc.  Thus, 
instead  of  regarding  the  multipath  phenomenon  as  a  nuisance  disturbance 
whose  effects  are  to  be  suppressed,  it  should  be  regarded  as  an  opportunity 
to  improve  system  performance. 

Two  bodies  of  work  in  the  literature  are  concerned  with  multipath 
receiver  design.  The  older  (see,  e.g.,  [22],  [27],  [23],  [l],  [4]) 
concentrates  on  the  explicit  diversity  structure  of  the  multipath;  its 
thrust  is  to  take  advantage  of  this  structure  bv  optimally  combining  the 
contributions  of  different  paths.  In  its  simplest  form,  this  approach 
assumes  resolvability  of  the  paths  (see  (4))  and  ignores  intetsymbol 
interference,  the  latter  assumption  often  being  justified  by  taking  the 
duration  of  a  transmitted  symbol  to  be  large  compared  with  the  duration 
of  the  multipath  profile. 

More  recently  [15],  [16],  [17],  [18],  equalization  techniques  that  were 
developed  for  data  transmission  over  telephone  lines  [14]  have  been  applied 
to  the  radio  multipath  problem.  Here,  receiver  design  concentrates  on 
reduction  of  the  effects  of  intersymbol  interference,  and  the  diversity¬ 
combining  properties  of  the  receiver  are  only  implicit.  This  approach 
appears  most  suitable  when  the  paths  are  not  resolvable  and  when  the 
symbol  duration  is  smaller  than  the  multipath  profile's  "spread." 

A  melding  of  the  two  approaches  would  appear  desirable  but  has  not 
been  accomplished.  Since  we  are  concerned  here  with  the  case  in  which 
resolvability  condition  (4)  is  satisfied,  we  shall  therefore  pursue  only 
the  former  diversity-oriented  approach,  as  modified  to  take  into  account 
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the  deleterious  effects  of  intersymbol  interference.  Instead  of  indulging 
in  general  and  complex  derivations,  however,  we  shall  present  results 
using  a  tutorial  "building  block"  approach,  using  intuitive  arguments 
that  are  justified  by  references  to  more  formal  developments  in  the 
literature. 


A.  The  Optimal  Single-Path  Receiver 

We  start  with  the  simple  case  in  which  the  channel  comprises  only 
one  path  —  K=1  in  (3).  We  assume  initially  that  the  path  strength  aQ 
and  delay  t  are  known  (to=0  fo r  simplicity),  but  that  the  carrier  phase 

* 

0o  is  unknown,  being  a  random  variable,  uniformly  distributed  over  [0,2ti). 

Since  the  absence  of  multipath  implies  the  absence  of  intersymbol 
interference  (a  point  we  return  to  later),  we  can  concentrate  on  the 
reception  of  a  single  symbol,  say  over  the  interval  0  <  t  <  T.  Knowledge 
of  this  interval  of  course  implies  some  sort  of  synchronization  procedure 
at  the  receiver,  a  question  discussed  below. 

Suppose  the  received  signal  r(t)  is  as  in  (2),  with  0  t  ^  T,  and 
where,  in  (3),  K=l,  to=0,  aQ  Is  known,  and  9o  is  random  as  described  above. 
The  transmitted  signal  s(t)  of  (1)  can  be  any  of  M  possible  waveforms 

jw  t 

si(t)  =  Re[oi(t)  e  ],  i  =  1,...,M.  (5) 


We  assume  that  the  transmitter  has  chosen  among  the  s^  at  random  with 


equal  probability  and  that  the  have  equal  energy: 


T 

(  s^(t)dx  =  £  fnr  all  i. 

Jo  1 


The  additive  noise  n(t)  of  (2)  is  for  simplicity  assumed  to  be  white 

and  Gaussian,  although  non-Gaussian  noise  is  also  common  in  urban  radio 
* 

Random  path  phases  are  assumed  throughout  this  paper,  since  these 
generally  change  too  rapidly  in  the  mobile  environment  to  make  use  of 
coherent  receiver  techniques. 
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commun  i  ca  t  ion . 

It  is  well  known  that  the  optimal  receiver  —  i.e.,  the  receiver 
that  decides  which  was  sent  with  minimum  probability  of  error  —  has 
the  form  depicted  in  Fig.  5  (see,  e.g.,  [36]).  As  shown  there,  r(t)  is 
passed  through  a  bank  of  M  filters,  "matched"  respectively  to  s^(t), 
i*  1,...,M,  i.e.,  having  impulse  responses  Sj(T-t),  0  £  t  <■  T  [29], 

The  filter  outputs  are  envelope  detected  and  th  >-  envelopes  sampled  at 
t*T  and  compared.  The  index  i  =  1,...,M  of  the  largest  sample  is  the 
receiver's  output. 

In  Fig.  5,  we  have  shown  the  outputs  of  the  envelope  detectors  when 
Sj(t)  is  the  transmitted  signal  and  when  the  received  noise  component 
n(t)  is  negligible,  assuming  that  the  s^  have  been  "well  chosen."  This 
latter  assumption  means  that  if  we  define  the  complex  correlation  function 


Yik(t)  "  j  0*(T)°k(T“t)dT,  i,k  =  1 . M,  | t I  £  T,  (7) 

then  [ 27] 

lyik(t)i  <<  2&  all  k  j*  i,  all  i  (8a) 

I  Yii  (t )  1  <<  2^  for  1 1 1  >  all  i,  (8b) 

where  W  is  the  bandwidth  shared  by  nil  s^,  and, optimally  but  not  necessarily, 

Y.,(0)  =  0,  all  k  Mi  all  i  (8c) 

Ik 


In  sketching  the  envelope  detector  outputs,  we  have  also  assumed  that 

TW  >>  1,  i.e.,  the  signals  are  of  the  so-called  spread-spectrum  type  [6]. 

None  of  the  foregoing  assumptions  about  the  structure  of  the  signal  set 
M 

{s^(t)}^=^  is  necessary  for  the  optimality  of  the  single-path  receiver  of 
Fig.  5  to  hold;  but  we  shall  invoke  them  when  discussing  multipath 
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Figure  5.  Optimal  noncoherent-phase  receiver  for  single-path  channel. 
(Equiprobable,  equl-energy  signals,  Gaussian  noise.) 
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receivers  later,  as  they  become  necessary  or  desirable. 

The  nolsefree  waveforms  sketched  in  Fig.  5  are  in  fact  given  by  [27] 


ei(t)  ~  ilv(t~T)l’  *  =  1‘*‘-*M*  0 


t  <  2T. 


(9) 


Conditions  (8a, b)  and  TW  >>  ]  assure  that  the  j^th  output  envelope  e^(t) 
consists  of  a  sharp  "mainlobe"  peak  surrounded  by  low-level  "sidelobes," 
while  all  other  outputs  have  only  low-level  sidelobes.  By  careful 
signal  selection,  (8a, b)  can  be  satisfied  with  the  maximum  sidelobe 
level  in  all  these  waveforms  at  a  factor  of  about  2//fw  down  from  the 
mainlobe.  Typically,  for  TW=100,  this  means  that  the  maximum  sidelobe  is 
about  17  dB  or  more  down  from  the  mainlobe. 

If  condition  (8c)  is  also  satisfied,  the  values  of  e^(T)  for  Jtjf j  are 
zero  at  the  sampling  instant  t=T,  so  that  —  in  the  absence  of  receiver 
noise  —  the  receiver  will  not  make  an  error.  For  moderate  amounts  of 
receiver  noise,  the  probability  that  the  jtth  output  exceeds  the  j^th  at  t=T 
for  some  Z^i  is  very  small. 

A  final  feature  of  Fig.  5  is  important.  There,  we  have  depicted 
the  output  waveforms  when  a  single  isolated  symbol  is  sent  during 
0  <  t  <  T.  If  another  symbol  —  say  s^(t)  —  were  sent  immediately 
afterward,  in  T  £  t  <  2T,  it  is  clear  that  the  response  to  it  would  occur 
over  the  interval  T  <  t  _<  3T.  The  mainlobe  peak  in  the  i_th  output  would 
be  centered  exactly  at  t-2T,  precisely  when  all  responses  from  the  first 
symbol  have  died  out.  Thus,  on  sampling  the  outputs  at  t=2T,  one  would 
be  able  to  make  a  decision  based  on  the  response  to  the  second  symbol 
alone,  whence  our  previous  statement  that  no  intersymbol  interference 


occurs  in  this  case. 


If  we  should  attempt  to  use  the  receiver  of  Fig.  5  when  many  paths 


are  present  (K  >  1  in  (3)),  we  would  expect  from  the  linearity  of  the 
medium  and  of  the  matched  filters  that  the  envelope  detector  output  wave¬ 
forms  will  look  something  like  those  in  Fig.  6.  Here,  we  have  shown  a 
four-path  situation  (K=4). 

The  £th  response  in  Fig.  6  is  the  envelope  of  the  superposition  of 
the  several  paths'  contributions,  and  —  when  noise  is  absent  —  can  be 
shown  from  (3)  and  (7)  to  be  of  the  form  [27, 


e»(t) 


A  1 


K-l  3\  * 

£  a  e  \  o  (t)o  (r+T+t  -t)di  , 
k«=0  J  0  J 


l  *  1,... ,M,  0  <  t  <  2T  +  A. 


(10) 


Under  resolvability  condition  (4),  the  mainlobe  peaks  in  the  jth  output 
ej(t)  are  distinct,  and  occur  as  shown  at  tQ  =0,  t  ^ ,  t£»  and  t^  ■  A. 

(The  quantity  t^  ^  -  tQ  =  A  is  called  the  multipath  spread;  it  is  by 
this  amount  that  the  waveforms  of  Fig.  6  spread  beyond  those  of  Fig.  5.) 

The  heights  of  these  peaks  are  proportional  to  the  path  strengths  a^. 

The  sidelobes,  both  of  ej (T)  and  of  the  other  outputs  (none  of  the  latter 
having  mainlobe  peaks),  are  mixtures  of  sidelobes  due  to  the  several  paths. 
We  stress  that  Fig.  6  is  drawn  for  the  isolated  transmission  of  a  single 
waveform  Sj(t),  0  _<  t  £  T. 

The  waveforms  of  Fig.  6  differ  from  those  of  Fig.  5  in  several 
important  respects: 

(a)  Strong  peaks  are  available  in  ej  (t)  at  multiple  times.  If  the 
decision  circuit  of  Fig.  5  knows  the  values  of  the  path  delays 
t  ,...,t„  .,  it  can  sample  the  contributions  of  all  paths  and 

O  X 


combine  them,  affording  the  receiver  the  advantages  of  diversity 
reception,  as  discussed  earlier. 

(b)  The  sidelobe  levels  of  all  outputs  is  increased,  since  (10) 
shows  the  addition  of  multipath  contributions. 

(c)  The  responses  to  the  symbol  sent  during  0  <_  t  £  T  now  extends 
beyond  t=2T,  thus  overlapping  with  the  responses  to  the  next  symbol, 
which  is  sent  during  T  <_  t  _<  2T.  That  is,  we  now  have  intersymbol 
interference,  caused  by  the  multipath. 

Effects  (b)  and  (c)  are  deleterious,  while  (a)  is  favorable.  As  we  shall 
see,  however,  the  benefits  of  (a)  usually  far  outweigh  the  deterioration 
caused  by  (b)  and  (c). 

For  the  time  being,  we  shall  ignore  the  effects  of  intersymbol 
interference,  and  inquire  into  the  structure  of  the  optimum  receiver 
for  reception  of  a  single  symbol  through  multipath,  assuming  first  that 

K_i 

the  path  delays  {t^Q  are  known.  However,  we  again  assume  random 
K— 1 

phases  >  independently  and  uniformly  distributed  over  [0,2tv);  we 

K-l 

also  assume  that  the  path  strengths  are  random,  perhaps  having 

different  distributions. 

Intuitively,  one  might  expect  under  these  conditions  that  the  optimal 
receiver  is  still  of  the  form  of  Fig.  5,  but  where  the  decision  circuit  now 
samples  each  of  its  inputs  at  multiple  times  T+t^,  k  =  0,...,K-1,  combines 
these  samples  for  each  input,  and  compares  the  resulting  combined  values; 
the  decision  would  be  the  index  of  the  largest  combined  value.  Indeed 
this  is  the  case,  at  least  when  (4)  and  (8b)  hold  so  that  the  pulses  in 
output  j  of  Fig.  6  are  distinct  [27].  However,  the  optimal  combining 
law  is  sometimes  complicated,  and  depends  on  the  statistics  of  the  path 


strengths . 


Suppose  that  the  sample  of  the  £_th  output  envelope  at  time  T+t^  is 
(In  the  absence  of  noise  x  =  a«  given  by  (10).)  Then,  if 

all  path  strengths  are  known,  the  optimal  combining  of  the  samples  is 
given  by  [27] 


^akXHk 

k=0  o 


where  I  is  a  Bessel  function  and  where  N  is  the  channel  noise  power 
o  o  v 

density.  If,  on  the  other  hand,  the  kth  path  strength  is  Rayleigh 

distributed  with  mean-square  strength  =  E[a^],  and  all  path  strengths 
* 

are  independent,  the  optimal  combining  law  is  [27],  [30] 


-  V  vik 


k=0  N  +0,  O 
o  k 

where  £is  the  common  energy  of  the  signals  s^,  given  by  (6).  More 
complicated  combining  laws  for  other  strength  distributions  are  given 
elsewhere  [3],  [27],  [30].  In  any  case,  a  decision  is  made  by  comparing 
the  Wj,  and  favoring  the  largest. 

Note  that  different  combining  laws  accentuate  the  various  samples 
in  different  ways.  In  (11),  for  example,  the  samples  are  approximately 
linearly  combined,  since  IQ(>)  increases  approximately  exponentially 
with  its  argument  for  large  argument;  but  samples  corresponding  to  paths 
with  larger  strengths  are  given  more  weight.  In  (12),  the  samples  are 
square-law  combined,  thus  accentuating  the  larger  samples;  but  all  samples 
from  paths  with  large  mean-square  strength  (ip^Q>>No)  are  weighted 
equally  while  samples  from  weaker  paths  are  suppressed.  The  essence  of 

* 

The  assumption  of  independent  path  strengths  is  a  variation  from  the 
multipath  model  described  in  the  previous  section,  and  to  that  extent 
the  resulting  receiver  is  only  quasi-optimal. 


optimal  combining  laws  is  the  relative  accentuation  of  more  credible 
data  and  the  relative  suppression  of  less  credible  data. 
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C .  The  Optimal  Multipath  Receiver:  Unknown  Delays 

K— 1 

As  indicated  in  Section  II,  the  path  delays  { ^  and  the  number 
of  paths  K  are  often  random  variables,  nor  known  a  priori.  In  order  to 
determine  the  optimal  receiver  for  this  situation,  we  simplify  somewhat 
from  the  path-delay  model  described  in  Section  II.  We  now  assume  that 
the  are  independently  chosen  from  a  single  common  probability  density 

distribution  p(t^),  0  t^  ^  A,  and  their  indices  subsequently  re-ordered 
in  order  of  increasing  t^.  This  model  violates  the  assumptions  in 
Section  II  in  two  respects.  First,  the  resolvability  condition  (4)  is 
not  always  met,  since  it  is  possible  that  two  paths  will  be  drawn  from 
the  distribution  in  such  a  way  that  |t^-t^|  <  However,  the  probability 

that  this  will  occur  is  small  if  there  arc  no  intervals  of  length  £  ^ 
in  which  substantial  probability  is  concentrated;  so  pCt^)  must  be 
"diffuse,"  without  high  peaks  and  with  WA  >>  1.  Second,  the  method  of 
generation  of  strings  discussed  in  Section  II  (see  Fig.  2)  incorporates 
dependences  among  i^'s  for  neighboring  i's,  which  implies  corresponding 
dependences  in  the  associated  delays  that  are  not  incorporated  in  the 
simplified  model  just  broached. 

These  variations  from  reality  are  not  substantial  for  purposes  of 
deriving  a  receiver  that  will  he  quas i-optimal  for  the  real  channel. 

The  density  distribution  p(t.  )  of  the  simplified  model  can  be  determined 
from  the  empirical  data  described  in  Section  II:  it  is  just  the  path- 

It  is  tempting  to  apply  subsequent  results  to  a  "channel  sounding" 
receiver,  in  which  pCt^)  becomes  an  a  posteriori  distribution.  But  such 
a  distribution  would  be  highly  peaked,  and  would  therefore  violate  this 
"dif fuseness"  condition.  We  shall  comment  further  on  this  point  later, 
when  discussing  channel-sounding  receivers. 
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occupancy  curve  exemplified  by  Fig.  4(a),  as  normalized  to  unit  area  by 
dividing  it  by  the  average  number  of  paths  E(K).  Were  the  simplified 
model  used  to  generate  tiie  T ^  strings  of  Fig.  2,  the  path-number 
distributions  would  be  Poisson  distributions  instead  of  the  somewhat 

•k 

narrower  distributions  exemplified  in  Fig.  4(b). 

In  deriving  the  quasi-optimal  receiver  for  unknown  delays,  we  assume 
that  path  strengths  a^  and  a^  are  independent  for  all  MU,  a  deviation 
from  the  reality  that  paths  whose  delays  are  not  greatly  different 
generally  have  correlated  strengths.  As  with  our  simplified  delay  model, 
we  follow  our  intuition  in  assuming  that  the  derived  receiver,  based  on 
the  simplified  strength  model,  will  be  close  to  optimal  in  the  real  world. 

With  these  mathematical  simplications,  and  assuming  that  (8b)  holds, 
the  optimal  receiver  structure  can  be  easily  derived  [27],  [4].  This 
receiver  computes  the  quantities 

a  rT+A 

W£  J  p(t-T)F[xe(t),t]dt,  1  «  1.....M,  (13) 

where  p(-)  is  the  path-delay  density  defined  above,  x„(-)  is  the  output 
envelope  of  the  £_th  filter,  and  F  [  - ,  -  ]  is  a  time-varying  nonlinear  function. 
The  w^'s  are  compared  and  a  decision  made  favoring  the  index  of  the 
largest  w  . 

The  nonlincarities  F[*,*]  depend  on  the  path-strength  statistics. 

If  the  path  strengths  are  known,  then  (cf.  (11)) 

2a(t-T)Xp (t) 

F[x4(t),t]  =  I  [ - jj - 1  (14) 

o 

where  a(t)  is  the  strength  of  a  path  .at  delay  L.  If  the  path  strengths 

*See  [33],  Fig.  5,  for  comparisons  of  these  Poisson  distributions  with 
the  actual  empirical  distributions. 


are  all  Rayleigh  distributed,  with  a  path  at  delay  t  having  mean-square- 
strength  <f/(t),  then  (cf.  12)) 
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<P(t-T)x  (t) 

F[x  (t),t]  =  exp[ - Tf]  (15) 

*  N  +i]/(t-T)0 

Expressions  for  F f - , * ]  for  other  path-strength  statistics  are  given  in 
[3],  14],  [27]. 

In  general,  F[x,t]  is  positive  and  monotone  increasing  in  x,  and  can 
therefore  be  written  as 


F[x,t]  =  F[ 0, t]  +  F[x,t], 


where  F[x,t]  is  positive  and  monotone  increasing  in  x.  On  substitution 
of  (16)  into  (13),  the  F[0,t]  term  leads  to  a  component  of  w£  that  is 
independent  of  1  and  can  therefore  be  neglected  in  the  comparison  of  the 
w£'s.  That  is,  the  receiver  need  only  compare  the  quantities 


T+  A 

w£=j  p(t-T)F[xe(t),t]dt,  i  =  1 . M. 


We  note  that  these  integrals  can  be  realized  by  convolution,  i.e.,  by 
generating  the  functions 


ve(t) 


p[A-(t-T)]F[x£(x),T]dT 


and  sampling  them  at  t  =  T+A.  Thus,  w£  can  be  realized  by  passing  x£(t) 
into  a  nonlinearity  F,  passing  the  output  of  F  into  a  filter  with  impulse 
response  h^(t)  =  p(A-t),  and  sampling  the  filter's  output  at  t  =  T+A. 

The  decision  circuit  of  Fig.  5  is  therefore  replaced  by  the  circuitry  of 


Fig.  7. 


Notice  that,  in  contradistinction  to  the  combining  laws  of  (11)  and 
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Figure  7.  Optimal  noncoherent  receiver  for  multipath  channel  with 
unknown  delays.  (Equiprobable,  equi-energy  signals; 
Gaussian  noise;  Poisson  path  delays,  independent  path 
strengths . ) 
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(12),  in  whic  the  samples  x^»  k  =  0,...,K-1,  are  combined  linearly 


or  quadratically ,  the  combining  law  of  (18)  involves  extreme  nonlinearities 
of  the  exponential  type.  For  example,  if  the  path  strengths  are  all 
Rayleigh  distributed  with  a  common  mean-square  value  i|-,  then  the 
nonlinearity  in  (18)  is  time-invariant,  and  of  the  form 


Flx^.t]  =  exp [ 


N  +!//< 
o 


-  1, 


(19) 


which  is  shown  in  Fig.  8.  When  passed  through  this  nonlinearity,  the 
1th  envelope  detector  output  of  Fig.  6  is  transformed  into  the  waveform 
of  Fig.  9.  (On  the  same  scale,  the  transformations  of  the  other  waveforms 
of  Fig.  6  are  negligibly  small.) 

The  waveform  of  Fig.  9  illustrates  a  certain  self-adaptivity  implicit 
in  (17).  A  priori,  path  delays  are  unknown.  However,  a  large  pulse  in 
one  of  the  x«(t)*s  at  some  instant  t  in  [T,T+A]  is  convincing  evidence 
that  a  path  exists  at  delay  t-T  (see  Fig.  6);  the  larger  the  pulse,  the 
more  convincing  is  the  evidence,  and  the  more  heavily  the  pulse  is 
emphasized  by  the  nonlinearity.  On  the  other  hand,  whenever  x,(t)  is 

/i 

small,  it  is  presumed  to  be  caused  by  noise  (or  sidelobes)  and  it  is 
strongly  suppressed  by  the  nonlinearity.  The  output  of  the  nonlinearity 
thus  presents  data  that  are  heavily  adjusted  by  a  posteriori  evidence  of 
the  existence  of  paths.  As  shown  by  (17),  these  data  are  further 
weighted  by  the  a  priori  knowledge  of  the  probabilities  of  path 
occurrences  implicit  in  the  function  p(-). 

We  stress  that  the  illustrations  in  Figs.  6  and  9  are  drawn  for 
relatively  large  average  signal-to-noise  ratio  ^£/No,  for  which  case  the 
signal  peaks  are  prominent  and  are  highly  emphasized  with  respect  to  the 
noise  by  the  operation  of  F.  On  the  other  hand,  when  is  smaller, 
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the  noise-suppressing  effect  of  the  nonlinearities  will  not  be  great; 
that  is,  the  various  outputs  F[x^(t),t],  i  =  1,...,M  —  even  the  j_th  — 
will  be  of  the  same  scale.  The  signal  peaks  in  the  J[th  output  will  then 
contribute  comparatively  little  to  the  ^th  integral  in  (17),  and  the 
receiver  will  be  prone  to  error.  One  can  see  that  known-pa th-delay 
receivers  (such  as  these  based  on  (11)  and  (12)),  by  only  having  to 
sample  the  envelope  detector  outputs  at  the  positions  of  signal  peaks 
and  not  the  noise  contributions  at  other  instants,  will  perform  better 
than  unknown-path-delay  receivers. 

D.  A  DPSK  Receiver 

As  previously  mentioned,  phase-coherent  techniques  have  been  avoided 
here  because  of  the  complexity  of  coherent  receivers  and  because  of  the  rapid 
time  variations  of  path  phases  with  vehicle  motion.  On  the  other  hand, 
differentially  coherent  techniques  clearly  will  show  promise  if  the  path 
phases  do  not  vary  appreciably  over  the  interval  during  which  two 
successive  signals  are  sent;  this  is  the  usual  case.  Although  the  form 
of  the  optimal  multipath  receiver  for  differentially  coherent  signalling 
has  not  yet  been  rigorously  derived,  it  is  strongly  suggested  by  that 
for  the  one-path  case. 

Let  a  signal  of  the  form  of  s(t)  of  (1)  be  differentially  phase- 

shift  keyed:  a  binary  0  is  sent  by  following  a  previous  transmission 

—  say  +s(t),  t  £  [0,T]  —  by  a  transmission  of  the  same  polarity  in  the 

succeeding  interval  —  +s(t-T),  t  G  [T,2T];  a  binary  1  is  sent  by 

changing  polarity  —  +s(t)  followed  by  +s(t-T).  In  the  single-path  case, 

the  received  waveform  in  (0,T]  is  (cf  (l)-(3),  with  K=l,  t  =0) 

o 

jo)  t  je 

r^(t)  =  +a}  Re[o(t)e  °  e  °]  +  n^(t),  0  _<  t  <  T,  (20) 


r 


l 
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where  n^(t)  is  a  noise  waveform.  If  d-0  or  1  is  the  data  symbol  sent, 
the  received  waveform  in  (T,2T]  is 

Jco  (t-T)  JO 

r2(t)  =  ±(-D  ao  Re[o(t-T)e  e  ]  +  n^t),  (21) 

T  <  t  _<  2T, 

where  n„(t)  is  a  noise  waveform  and  we  have  assumed  that  a  and  0  have 
2  o  o 

not  changed  from  their  values  during  (0,T],  The  optimal  receiver  in  this 
case  is  well  known  [24]  to  have  the  form  of  Fig.  10.  (The  lowpass  filter 
in  Fig.  10  serves  only  to  eliminate  the  double-frequency  terms  generated 
by  the  multiplier.) 

By  performing  the  operations  shown  in  Fig.  10  on  the  signal  components 
of  (20)  and  (21)  —  i.e.,  neglecting  the  noise  components  —  one  can 
easily  show  that  the  output  of  the  lowpass  filter  has  the  form 

y(t)  =  \  (-l)V|y(t-2T)|2  (22) 

where 


and  y(t)  -  0  elsewhere.  On  comparing  (22),  (23)  with  (7),  (9),  we  see 
that  y(t)  is  twice  the  squared  envelope  of  the  response  of  the  matched 
filter  to  the  waveform  s(t),  delayed  by  T  seconds  and  keyed  by  (-1)^. 
Thus,  in  the  absence  of  noise,  the  decision  circuit  will  output  exactly 
the  input  digit;  when  noise  is  present,  errors  will  of  course  occur. 

We  now  conjecture  that  for  tiie  case  of  multipath  obeying  the 
simplified  model  of  Section  C  above  —  i.e.,  Poisson-distributed, 
resolvable  path  delays  and  independent  path  strengths  —  the  optimal 
DPSK  receiver  bears  the  same  relationship  to  Fig.  10  as  the  receiver  of 
Fig.  7  bears  to  that  of  Fig.  5.  More  precisely,  we  conjecture  that  the 


•f 


1 


t  \ 


optimal  receiver  has  the  form  of  Fig.  11.  There,  the  nonlinearity 
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F(y,t)  =  F(/[yj,t)  sgn  y  (24) 

is  a  bipolar  version  of  F,  adjusted  for  the  fact  that  the  nonlinearity's 

input  is  related  to  the  square  of  the  matched  filter  output  envelope, 

rather  than  —  as  in  Fig.  7  --  the  envelope  itself.  A  graph  of  F  for 

the  F  of  (19)  is  shown  in  Fig.  12;  compare  this  to  Fig.  8. 

The  samples  taken  every  T  sec  by  the  decision  circuit  in  Fig.  11 

are  timed  to  capture  the  extrema  of  the  output  of  the  path  integrating 

filter  hp(-)-  Figure  13  shows  some  appropriate  waveforms  illustrating 

this  point.  In  the  absence  of  noise,  the  output  y(t)  of  the  product 

detector  of  Fig.  10,  and  therefore  the  input  to  the  nonlinearity  in 

* 

Fig.  11,  is  approximately  of  the  form  (cf.  (22)) 


y(t)  =  |  (~l)d  £  a^jy(t-2T-tk)|Z,  2T  <  t  <  2T+A,  (25) 

where  A  is  —  as  before  —  the  multipath  spread.  A  sequence  of  such 
noiseless  y(t)'s  Is  shown  in  Fig.  13(c),  for  the  input  symbol  sequence 
10110...  The  corresponding  outputs  of  F  and  h^(-)  are  shown  in  Figs. 
13(d)  and  13(e),  assuming  that 


h  (t) 
P 


0  £  t  <  A 
elsewhere  . 


(26) 


In  (25),  we  have  assumed  that  the  sidelobes  of  y(t)  of  (23)  satisfy  a 
condition  of  the  form  of  (8a),  so  that  "interpath  interference"  is 
negligible;  i.e.,  the  sidelobes  due  to  path  l  arc  small  at  the  peak  due 
to  path  k,  i^k,  i  =  0,...,K-1.  We  have  also  neglected  intersymbol  inter¬ 
ference,  i.e.,  assumed  that  the  sidelobes  due  to  all  paths  in  one  signalling 
interval  are  negligible  insofar  as  they  extend  into  adjacent  intervals. 

We  return  to  the  questions  of  interpath  and  intersymbol  interference 
below. 


DECISION 

CIRCUIT 


Figure  11.  Conjectured  optimal  DPSK  receiver  for  multipath  channel 

with  unknown  delays. (Equiprobable  symbols;  Gaussian  noise; 
Poisson  path  delays;  independent  path  strengths.) 


Figure  12.  The  nonlinearity  of  Fig.  11  for  Rayleigh-distributed 
path  strengths  with  a  common  mean-square  value  ip. 
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The  output  of  hp(-)  is  sampled  at  instants  uT+A,  n  =  2,3,...,  and  leads 
to  the  receiver  outputs  shown;  those  are  delayed  by  T+A  seconds  from  the 
corresponding  input  symbols.  Of  course,  when  noise  is  present,  some  of 
the  output  symbols  will  differ  from  the  associated  input  symbols. 

E.  Channel-Sounding  Receivers 

Up  until  this  point,  our  discussion  of  receiver  design  has  been 
based  upon  a  priori  statistical  knowledge  of  the  channel.  In  many 
situations,  measurements  can  be  made  of  channel  characteristics  by  use 
of  sounding  signals  that  enable  the  derivation  of  a  posteriori  statistics. 
Such  sounding  signals  might  be  special  signals  used  for  sounding  only, 
or  might  be  the  data  signals  themselves;  in  the  latter  case,  sounding 
and  data  transmission  occur  simultaneously. 

One's  first  impulse  is  simply  to  use  the  receiver  structures 
discussed  above  for  the  case  of  unknown  delays,  but  to  base  the 
characteristics  of  F,  F  and  ( - )  in  Figs.  7  and  11  on  a  posteriori  rather 
than  a  priori  statistics.  This  would  in  fact  be  appropriate  with  regard 
to  the  path-strength  distributions  and  the  nonlinearities  F  and  F  they 
determine.  However,  as  indicated  in  the  footnote  on  page  23,  as  soon  as 
the  sounding  signals  enable  very  accurate  estimation  of  the  path  delays, 
the  diffuseness  condition  on  the  path-arival  distribution  p(t^)  no 
longer  holds,  and  the  resulting  filter  operation  h  (•)  that  is  based  on 
this  distribution  is  no  longer  even  quasi-optimal .  Although  expression^ 
for  optimal  receivers  using  nondiffuse  a  posteriori  delay  distributions 
can  be  derived,  they  and  the  resulting  receivers  become  inordinately  and 
unnecessarily  complicated. 

An  alternative  approach  is  usually  used.  In  this  approach,  it  is 
assumed  that  sounding  results  in  extremely  accurate  estimates  of  the  path 
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variables  (t^),  (B^),  or  at  lensc  °f  the  delays  and  strengths, 

if  not  the  phases.  These  estimates  are  then  assumed  to  be  exact,  and 
used  as  parameters  in  a  receiver  that  is  based  upon  exact  knowledge  of 
the  associated  variables. 

In  many  cases,  estimates  of  (0  }  are  not  deemed  worth  making,  either 
because  of  the  complexity  of  the  resulting  receiver  or  —  especially  in 
the  mobile  receiver  context  —  because  these1  phase  shifts  change  too 
rapidly  to  be  tracked  and  used  effectively.  Therefore,  we  discuss  iere 
only  noncoherent  channel-sounding  receivers  that  make  use  only  of  path- 
strength  and  path-delay  estimates. 

The  "optimal"  M'ary  strength/delay-cstimatlng  receiver  that  follows 
the  philosophy  just  described  is  clearly  based  on  (11).  The  quantities 


K-l  2a.  x„, 

wi  =  E  ioge  v-f**) 

*  k=0  N 


(27) 


must  be  calculated,  where  a^  is  the  estimate  of  the  strength  of  the  kth 
path  and  is  a  sample  of  the  J.th  matched  filter  output  envelope  at 

/\  a 

t  =  T+t^,  being  the  estimate  of  the  kth  path's  delay.  The  index 
( i =  1,...,M)  for  which  w^  is  maximum  is  the  receiver's  digital  output. 
If  the  paths  are  strong  enough  with  respect  to  the  noise  to  be  measured 
accurately  —  which  is  our  assumption  —  then  (27)  can  be  approximated 
by 

K-l 


w  =  —  y; 

l  N  _  M 


o  k=0 


akX£k 


(28) 


since  log  ^  I^(x)  s  x  f°r  lar8c  x-  This  approximation  is  the  optimal 
linear  diversity  combiner  of  Brennan  [2],  and  we  shall  use  it  henceforth. 

The  linear  combiner  of  (28)  can  be;  realized  through  use  of  a 
transversal  filter,  as  shown  In  Fig.  14.  This  filter  incorporates  a  delay 


1 


re  14 .  A  transversal  filter  used  in  realizing  (28) 
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line  A  seconds  long,  which  is  tapped  at  least  every  1/W  seconds,  for  a 
minimum  of  WA  taps.  The  input  to  the  transversal  filter  is  the  output 
envelope  of  a  matched  filter,  say  x^(t)  of  Fig.  6.  The  output  of  the 
transversal  filter  is  a  weighted  sum  of  certain  tap  outputs,  the  taps 
that  are  included  in  the  sum  depending  on  the  path  delay  estimates  t^. 

The  estimates  t^  (k=0, . . . ,K-1)  are  used  to  turn  on  the  amplifiers 
of  those  taps  having  the  delays  (measured  from  the  right-hand  end  of  the 
line 'i  that  most  closely  approximate  the  t^'s;  i.e.,  K  of  the  tap  amplifiers 
are  activated.  The  gains  of  the  activated  tap  amplifiers  are  then  set 
to  be  proportional  to  the  associated  strength  estimates  a^.  Amplifier 
gains  are  shown  in  Fig.  14  for  the  four-path  response  assumed. 

To  explain  how  the  transversal  filter  works,  we  have  shown  in  Fig.  14 
the  voltage  profile  along  the  delay  line  that  would  occur  if  the  delay 
line's  input  were  the  jth  matched  filter  output  envelope  x^ (t)  of  Fig.  6; 
this  profile  moves  to  the  right  with  time  and  is  shown  at  the  instant 
t  -  T+A.  At  this  instant,  the  peaks  in  the  profile  lie  at  delays 
t  ,t.,...,t  .,  as  measured  from  the  line's  right-hand  end.  If  the  delay 
estimates  tj.  are  accurate,  the  activated  taps  will  sample  the  profile 
close  to  these  peaks,  so  at  time  t  =  T+A  the  transversal  filter's  output 
will  be  approximately 

K-l  ^  K-l 

23  a  i  (T+t  )  *  £  a.x  .,  (29) 

k=0  J  k=0  1 

which  is  proportional  to  w^,  p-  =  j ,  of  (28). 

Actually,  the  transversal  filter  is  itself  a  matched  filter  of  sorts. 
Note  that  the  lowpass  equivalent  of  the  channel's  impulse  response,  i.e., 
that  relating  o(t)  and  p(t)  of  (1)  and  (3),  is 
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K.-1 


JO,. 


1.(0  =  £  a  e  «(t-t  ). 


k=0 


(30) 


On  the  other  hand,  the  impulse  response  of  the  transversal  filter  can 
clearly  be  written  as 


K-l 

=  E  a  6(t-A+t.) 
“  k=0  k  k 


K-l 

=  E  V(A"t"\)  on 

k=0 


The  transversal  filter  is  therefore  matched  to  an  estimate  of  the  magnitude 

* 

of  the  lowpass-equivalent  channel  impulse  response. 

As  Xj(t)  moves  to  the  right  in  Fig.  14,  the  transversal  filter  con¬ 
volves  Xj(t)  with  ht^(t),  producing  an  output  like  that  shown  in  Fig.  15. 
The  main  peak  of  this  output  occurs  when  t  =  T+A,  i.e.,  when  Xj(t) 
is  aligned  as  in  Fig.  14,  with  its  peaks  all  located  at  activated  taps. 

This  main  peak  is  proportional  in  height  to  (29).  The  output  also  has 

minor  side  peaks  proportional  to  a,x.  (k=0, 1 , 2 , 3 ;r^k)  and  sidelobes 

K  j  r 

accumulated  from  convolution  of  h  -(t)  with  the  sidelobes  of  x.(t). 

tf  j 

It  is  clear  that  we  need  to  sample  only  the  main  peak,  to  obtain  (29). 

A  complete  M'ary  channel-estimating  receiver  can  therefore  be 
depicted  as  in  Fig.  16.  Each  envelope  detector  output  x^(t),  l  =  1,...,M, 
enters  a  transversal  filter  of  the  form  of  Fig.  14,  the  parameters  of 
which  are  driven  by  the  estimates  {t^},  (a^)  available  from  a  sounding 
receiver.  At  time  t  =  T+A, the  transversal  filters'  outputs  are  sampled  and 
compared,  and  a  decision  made  favoring  the  index  of  the  largest.  Such 


As  discussed  later  in  Section  IV,  if  phase-coherent  techniques  were  used, 
the  transversal  filter  would  become  a  bandpass  filter  matched  to  an  estimate 
of  the  channel  impulse  response  itself.  The  cascade  of  the  receiver's 
signal  matched  filter  and  the  transversal  matched  filter  would  then  form  a 
receiver  whose  filter  is  matched  to  an  estimate  of  the  actual  received 
signal,  as  modified  by  the  channel. 
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Figure  15.  Illustration  of  the  operation  of  the  transversal  filter  of 
Fig.  14:  (a)  Filter  input;  (b)  Filter  impulse  response; 

(c)  Filter  output,  obtained  by  convolving  (a)  and  (b) . 


»  _  j  *'W..  . 


M ttWM* 


44 


a  receiver  lias  been  called  a  RAKK  receiver  [23]  because  of  the  tooth-like 
structure  of  the  taps  on  rue  transversal  filter's  delay  line. 

The  sounding  receiver  of  Fig.  lb  can  itself  be  structured  as  a 
cascade  of  a  matched  filter,  envelope  detector  and  tapped  delay  line, 
where  the  matched  filter  is  matched  to  a  known  waveform,  perhaps  containing 
more  energy  and  lasting  longer  than  the  data  signals  s  .  Prior  to  data 
transmission  (and  during  it,  if  the  channel  changes  rapidly  enough), 
the  sounding  receiver  "listens"  for  its  signal.  When  the  signal  is 
received,  the  matched- f i 1  ter  output  envelope  will  have  multiple  peaks, 
like  the  jth  waveform  of  Fig.  h.  This  waveform  will  progress  down  a 
tapped  delay  line  like  that  in  Fig.  14.  A  threshold  is  triggered  when 
the  first  peak  readies  the  line's  right  end,  at  which  moment  the  voltages 
at  each  tap  are  sampled  and  held  if  they  exceed  another  threshold,  or 
set  to  zero  i f  they  do  not.  The  frozen  tap  voltages,  which  are  propor¬ 
tional  to  tile  path  strengths  at  the  tap's  delays,  are  then  used  to  set 
the  gains  of  the  associated  taps  in  the  data  receiver's  transversal 
filters . 

A  self-adaptive  version  of  this  receiver  uses  the  data  signals 

themselves  for  sounding.  At  the  time  of  a  decision  favoring,,  say,  signal 

j,  the  tap  voltages  of  the  jth  transversal  filter  are  frozen,  just  as  in 

the  sounding  receiver  discussed  above;  those  voltages  are  then  used  to 

* 

set  the  corresponding  tap  gains  of  all  the  transversal  filters.  This 
self-adaptive  version  can  be  started  by  using  a  threshold  triggering 
device  of  the  type  used  in  the  separate  sounding  receiver. 

We  note  that  the  sounding,  mechanism  also  performs  a  receiver 

synchronizing  function,  which  --  ns;  previously  mentioned  --  wc  have  thus 
k 

More  elaborately,  the  nth  tap  gain  can  he  set  on  the  basis  of  a  weighted 
average  of  tin-  nt!i  lap  voltages  held  over  a  given  number  of  past  decisions. 
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far  ignored.  Triggering  of  the  threshold  on  the  rightmost  tap,  in 
either  the  separate-sounding  or  self-adaptive  realization,  starts  a 
clock  that  subsequently  supplies  properly  timed  sampling  pulses  to  the 
decision  circuit. 

Finally,  we  note  that  a  RAKE-1  ike  DPSK  receiver  can  be  structured 
along  the  same  principles.  This  is  shown  in  Fig.  17.  Notice  that,  since 
the  input  in  Fig.  17  is  akin  to  the  square  of  the  inputs  of  Fig.  16,  we 
have  inserted  a  bipolar  square-root  operation  before  the  transversal 
filter  in  the  DPSK  realization. 

F.  Intersymbol  Interference 

Heretofore,  our  discussion  has  ignored  the  effects  of  intersyrabol 
interference  by  concentrating  on  isolated  single  transmissions 
(see,  e.g..  Fig.  6).  Even  in  the  DPSK  case,  where  information  is  sent 
via  the  agreement  or  disagreement  of  the  polarities  of  two  successive 
transmissions,  these  transmissions  were  assumed  to  be  sufficiently  spaced 
so  that  the  time  dispersion  caused  by  the  multipath  channel  causes  little 
or  no  performance  degradation  (see  Fig.  13). 

Suppose  now  that  each  transmission  of  a  choice  from  the  signal  set 
requires  T  seconds,  and  that  successive  transmissions  occur  at  intervals 
of  length  T^.  As  indicated  by  (10),  the  outputs  of  the  matched  filters 
in  the  receivers  discussed  above  will  last  2T+A  seconds  in  response  to 
each  transmission  of  length  T,  raising  the  possibility  of  intersymbol 
interference.  In  order  to  help  visualize  intersymbol  interference  in 
the  spread-spectrum  (TW>>1)  case  of  importance  here.  Fig.  18  shows 
sequences  of  output  envelopes  of  a  single  matched  filter  in  response  to 
a  periodic  input  of  the  signal  to  which  it  is  matched.  This  diagram 
illustrates  the  inter-relationships  among  T,  Tfi  and  A. 
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Figure  17.  A  channel-sounding  DPSK  receiver. 
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Figure  18.  Illustration  of  intersymbol  interference  (a)  No 

intersymbol  interference:  T  >_  T+A.  (b)  Moderate 
intersymbol  interference:  T  =  T,  A  =  2T. 
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In  Fig.  18(a),  we  have  shown  a  case  in  which  there  is  no  intersymbol 
interference,  a  situation  requiring  that  Tg  ^  T+A.  For  reference,  we  have 
also  shown  a  RAKE  delay  line  and  combiner,  of  length  A;  this  is  turned 
end  for  end  compared  to  that  of  Fig.  14.  The  output  envelopes  of  the 
response  of  a  matched  filter  to  three  successive  transmissions  are  shown, 
which  should  be  visualized  as  entering  the  delay  line  from  the  right  and 

i k 

moving  with  time  to  the  left.  This  output  sequence  is  shown  "frozen"  at 
the  instant  at  which  the  multipath  pulses  of  the  central  (nth)  transmission 
are  perfectly  aligned  with  the  delay  line  and  ready  to  be  sampled. 

Notice  that,  at  this  instant,  neither  the  (n-l)st  nor  the  (n+l)st  response 
is  in  the  delay  line  to  interfere  with  the  nth. 

In  Fig.  18(b)  a'  case  of  moderate  intersymbol  interference  is  shown. 
Here  =  T  and  A  =  2T,  so  Tg  =  ^  (T+A).  A  sequence  of  output  envelopes 
in  response  to  seven  successive  transmissions  is  shown,  frozen  at  the 
instant  when  the  multipath  pulses  of  the  central  (nth)  response  are 
aligned  with  the  delay  line.  However,  in  this  case,  parts  of  two 
predecessor  and  two  successor  responses  (n-2 , n-1 , n+1 , n+2)  are  also  in 
the  delay  line,  causing  intersymbol  interference.  In  general,  a  total 
of  2[t/a]  predecessor  and  successor  symbol  will  interference  with  each 
symbol  as  it  is  sampled,  whore  J” x |  is  the  smallest  integer  larger  than 
X. 


k 

Of  course,  in  practice  these  responses  would  have  been  superposed  by  the 
matched  filter  prior  to  envelope  detection,  and  the  actual  output  envelope 
would  be  a  nonlinear  combination  of  the  three  responses. 

kk 

In  the  DPSK  case  shown  in  Fig.  13,  T  =  T  and  A/T  <  1,  so  there  is 
intersymbol  interference  by  only  one  successor  and  one  predecessor. 

The  interference  is  due  only  to  the  sidelobes  of  the  multipath 

pulses  in  neighboring  transmission  intervals,  rather  than  to  the  pulses 

themselves  as  in  Fig.  17(b),  and  the  interference  can  be  deemed  negligible. 
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In  order  completely  to  avoid  Intersymbol  interference,  the  condition 
Tg  >  T+A  must  be  satisfied.  Although  we  have  some  flexibility  in 
decreasing  T  (subject  to  our  requirement  TW  >>  1  and  to  limitations 
imposed  on  the  bandwidth  W) ,  A  is  fixed  by  the  channel.  Thus,  for  a 
binary  system,  we  cannot  completely  avoid  intersymbol  interference  when 
the  data  rate  is  greater  than  roughly  1/A  bps.  For  example,  for 
binary  transmission  through  urban  multipath  (A  =  5  ps,  effectively),  any 
transmission  rate  greater  than  about  200  Kbps  will  result  in  intersymbol 
interference,  even  in  a  spread-spectrum  system. 

Of  course,  it  is  well  known  that  the  effects  of  intersymbol 
interference  can  be  ameliorated,  so  higher  data  rates  can  be  achieved 
without  undue  deterioration  of  receiver  performance.  One  approach, 
followed  by  Monsen  (15],  [16],  [17],  is  based  on  classical  equalization 
techniques  developed  for  reduction  of  intorsymbol  interference  on  baseband 
landlines;  it  is  particularly  applicable  when  TW  =  1  and  A  >>  T.  In  our 
case,  when  TW  >>  1  and  A/T  is  moderate,  we  follow  a  different  approach, 
based  on  the  RAKE  receiver. 

The  basis  of  the  RAKE  approach  is  to  recognize  that,  while  multipath 
pulses  from  predecessor  and  successor  symbols  are  on  the  RAKE  delay  line 
at  the  decision  instant  for  the  present  symbol, as  shown  in  Fig.  18(b), 
it  is  unlikely  that  they  will  appear  at  taps  which  are  activated.  Recall 
that  the  delay  line  has  taps  every  1/W  seconds.  Since  TW  >>  1,'  and  we 
are  now  assuming  that  T/A  >  1,  we  will  have  WA  >>  1.  There  will  thus  be 
large  number,  WA,  of  taps.  For  example,  if  TW  =  100  and  T/A  =  5,  then 
there  are  W  =  500  taps  on  the  delay  line.  On  the  other  hand,  only  those 
taps  at  which  pulses  are  expected  are  activated,  so  even  if  there  are 
as  many  as  20  paths  in  the  example  above,  only  20/500  =  2.5%  of  the  taps 
will  be  activated.  These  activated  taps  will  be  perfectly  aligned  to 
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sample  the  multipath  pulses  of  the  central  (nth)  response  in  Fig.  18(b); 
but  it  is  extremely  unlikely  that  the  pulses  of  adjacent  responses  on  the 
line  will  also  be  aligned  with  the  activated  taps. 

Yet  another  insight  into  the  capability  of  the  RAKE  transversal 
filter  to  suppress  intersymbol  is  given  by  looking  at  the  filter's  output. 
Although  the  total  input  to  the  delay  line  is  really  a  nonlinear  combination 
of  the  pulse  trains  shown  in  Fig.  18(b),  for  simplicity  one  can  visualize 
the  result  of  each  pulse  train  sweeping  to  the  left  through  the  line  and 
being  convolved  with  the  RAKE  filter's  impulse  response,  as  depicted 
in  Fig.  15.  A  sequence  of  seven  such  resulting  convolutions,  when 
properly  combined,  would  look  something  like  Fig.  19.  Each  of  the  major 
peaks  there  corresponds  to  the  exact  alignment  of  one  of  the  pulse  trains 
with  the  activated  taps  on  the  delay  line,  and  it  is  these  peaks  that  are 
sampled.  The  pedestal  upon  which  the  major  peaks  sit  is  composed  of  minor 
peaks,  as  in  Fig.  15,  and  a  general  sidelobe  "hash"  level.  We  call  this 
pedastal  "multipath-induced  interference." 

The  ability  of  the  RAKE  receiver  to  concentrate  on  the  part  of  the 
matched  filter  response  due  to  the  current  symbol  is  peculiar  to  its 
tapped-delay-1 ine  structure.  The  integrating  receivers  of  Figs.  7  and  11, 
by  integrating  over  the  interleaved  responses  to  several  symbol s, cannot 
suppress  intersymbol  interference  as  RAKE  does. 

In  summary,  it  appears  that  data  rates  much  greater  than  1/A  bps  can 
be  achieved  with  RAKE  receivers  without  undue  deterioration  of  performance. 
We  investigate  this  possibility  both  analytically  and  by  simulation  in 
the  next  several  sections. 

We  close  this  section  by  noting  one  complication  that  arises  in 
realization  of  a  RAKE  receiver  when  T/A  >  1.  Recall  that  one  alternative 
to  estimation  of  the  path  delays  and  strengths  was  self-adaptive,  using 
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the  data  signals  themselves  as  sounding  signals, and  using  the  multipath 
profiles  that  are  stretched  out  along  the  HAKE  delay  lines 
to  aid  in  path-parameter  estimation.  As  shown  in  Fig.  18(b),  however, 
when  the  present  (nth)  signal's  multipath  profile  is  aligned  with  the 
delay  line,  adjacent  signals  cause  spurious  pulses  to  appear  on  the  line. 
These  latter  must  be  ignored  in  order  for  the  RAKE  receiver  to  avoid 
intersymbol  interference,  but  —  in  the  adaptive  approach  to  path- 
parameler^estimation  —  they  would  tend  to  be  identified  as  true  paths  and 
cause  spurious  activation  of  taps.  Thus,  when  T/A  >  1,  the  self-adaptive 
approach  to  estimation  cannot  be  used,  and  separate,  isolated  sounding 
signals  and  an  associated  sounding  receiver  are  necessary. 
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IV .  Analysis  of  Mul tipath-Indueed  Interference 

We  have  just  concluded  by  a  number  of  heuristic  arguments  that, 
even  if  T  <<  T+A,  a  RAKE  receiver  will  work  satisfactorily,  despite 

•k 

multipath- induced  interpath  and  intersymbol  interference.  In  the  present 
section  we  present  an  analysis  that  bolsters  our  confidence  in  this 
conclus ion . 


A.  A  Simplified  Model 

For  simplicity,  we  consider  a  coherent  receiver,  since  the  noncoherent 
receivers  discussed  in  Section  III  lead  to  complicated  nonlinear  analysis. 
More  specifically,  we  analyze  a  binary,  coherent  PSK  system  whose  receiver 
filters  are  shown  In  Fig.  20.  Here,  a  bandpass  matched  filter  is  followed 
by  a  bandpass  transversal  filter,  rather  than  (as  previously)  by  a  cascade 
of  an  envelope  or  DPSK  detector  and  a  lowpass  transversal  filter.  As  in 
the  case  of  the  lowpass  transversal  filter  of  Fig.  14,  the  tap  gains  in 
Fig.  20  are  set  by  estimates  of  path  strengths.  Now,  however,  these  gains 
not  only  weight  the  tap  outputs  by  estimates  a^  of  the  path  strengths, 
but  also  try  to  correct  for  the  paths’  phase  shifts  0^  by  phase  shifting 
the  tap  outputs  in  the  opposite  directions,  by  amounts  -0^.  We  have 
written  the  non-zero  gains  in  Fig.  20  in  complex  form: 

-A 


-*  -  j  k 

\  \  e 


(32) 


The  lowpass-equi valent  impulse  response  of  the  transversal  filter  is 
clearly  of  the  form  (cf.  (31)) 


htf(t) 


>  * 


K-l 

£  a.  6 (A-t-t, ) 
k=0  k  K 


(33) 


See  footnote,  page  33. 
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We  therefore  see  that  ht^.(t)  would  he  matched  to  the  channel's  impulse 
response,  given  by  (30),  if  the  estimates  were  perfect,  i.e.,  if  ak  -  a^, 

=  °k’  and  *k  =  tk’  a11  k* 


We  continue  in  lowpass-equivalent  form  in  the  frequency  domain.  Let 

o(t),  0  <  t  <  T,  be  the  lowpass  equivalent  of  a  signal  transmitted  through 

the  channel  (of.  (1))  and  let  the  Fourier  transform  of  o(t)  be  S(f). 

Disregarding  channel  noise  for  the  moment,  the  channel  output  then  has 

spectrum  S(f)  H  (f),  where  H  (t)  is  the  Fourier  transform  of  h  (t)  of 
mm  m 


(30),  i.e.. 


K-l 

Hm(f)  =  £  "k 
k=0  k 


-j  2nf  t. 


(34) 


j°k 

■where  =  a^  e  .  The  matched  filter  in  Fig.  20  has  lowpass- 
equivalent  impulse  response  o(T-t)  and  transfer  function  S*(f)e 
so  its  output  has  spectrum  |s(f)|2  ( f ) e  j f T ^  Tjle  transversa^  filter 

has  transfer  function 


K-l  -j2irf  (A-t.  ) 

=  E  «k  e  k  (35) 

k=0  K 


( the  Fourier  transform  of  (33)),  whence  its  output  has  spectrum 


W(f)  =  { S  (f )  |  2  Hm(f)  Htf(f) 


=  |S(f) 


K-l  K-l 

E 

k=0  1=0 


C‘ka£ 


e-j2nf(T+A) 


(36) 


We  assume  for  simplicity  that  the  channel-sounding  procedure  has  led 
to  very  good  estimates,  so  we  set  a^  =  a^,  0^  =  (hence,  =  a^),  and 
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t,  “  t,  .  We  also  shift  the  time  origin  to  the  right  by  T+A  seconds 
k  k 

to  eliminate  the  factor  in  (36).  Then  (36)  becomes 


-  K-l  K-l 

w(o  =  i s(f ) t  52  Z  iv 

k=o  e=o 


* 


j2nf  (tk-tj ,) 


(37) 


The  corresponding  time  waveform  is  the  inverse  Fourier  transform  of  W(f), 
i.e. , 

K-l  K-l  * 

w(t)  =  X,  jL,  Y(t-t  +t  ),  1 1 J  <_  T+A  ,  (38) 

k=0  £=0  k  k  1 

where  y(t)  is  the  inverse  Fourier  transform  of  |s(f)|  ,  and  is  given 
by  (23). 

For  the  case  in  which  a(t)  is  sent  only  once  through  the  channel, 

(38)  is  the  lowpass-equivalent  output  of  the  transversal  filter,  shifted 

so  that  its  peak  lies  at  t=0.  (The  true  bandpass  output  is  ~  Re[w(t)e^27Tft] 

and  the  output  envelope,  —  |w(t)j,  has  a  shape  like  that  depicted 

in  Fig.  15(c),  where  the  central  peak  is  now  centered  at  t=0.)  If  we 

OC 

now  suppose  that  a  long  phase-shift-keyed  sequence,  52  d  o(t-nT), 

n--°»  n 

dn  =  +1,  is  sent  through  the  channel,  the  corresponding  lowpass-equivalent 
output  of  the  transversal  filter  will  be 

CP 

v (t 1  —  52  d  w(t-nT) .  (39) 

n=-°° 

The  output  envelope,  ~  jv(t)[,  will  now  look  like  Fig.  19,  where  the 
result  of  only  seven  transmissions  is  shown.  The  central  peak  in  Fig.  19 
would  correspond  to  the  n=0  term  in  (39),  now  centered  on  t=0. 


*  As  seems  obvious,  and  is  shown  in  Appendix  C,  if  the 
signal's  energy  is  much  greater  than  the  data  signals 
a  reasonable  assumption. 


channel- sounding 
energy,  this  is 


•  \ 
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The  analysis  in  this  section  is  concerned  with  the  degradation  in 
performance  caused  by  the  multipath— induced  interference.  Ideally,  this 
degradation  should  be  measured  in  terms  of  error  probability.  However, 
the  multipath  statistics  that  govern  the  error  probability  are  so 
complicated  as  to  make  this  ideal  goal  unattainable.  Instead,  we  shall 
perform  an  "equivalent  noise"  analysis. 

We  concentrate  on  the  central  peak  in  v(t),  i.e.,  on 


v (0)  =  £  dn  w(-nT) .  (41) 

n=-oo 

Noting  from  (38)  that  w(t)  is  nonzero  only  in  |tj  T+A,  we  see  that  (41) 
only  has  2N+1  nonzero  terms,  where  N  is  the  largest  integer  less  than 
(T+A)/T.  Thus,  using  (38),  (41)  becomes 


v(0)  =  £ 


n=-N 


K-l  K-l 

E  E 

k=0  0 


Vl  Y(nT-V^). 


(42) 


This  is  the  signal  component  of  the  transversal  filter  output  at  t=0. 
The  channel-noise  component  —  call  it  n(0)  —  adds  to  this,  leading  to 
a  total  output 


z (0)  =  v(0)  +  n(0) . 


(43) 


Our  approach  to  evaluating  the  effect  of  multipath-induced  interference 
is  to  determine  how  much  its  presence  increases  the  variance  of  2(0) 
above  the  value  contributed  by  n(0),  i.e.,  how  much  "equivalent  channel 
noise"  is  added  by  intersymbol  interference. 

In  order  to  proceed,  we  make  a  number  of  simplifying  assumptions 


about  the  multipath  statistics  that  are  akin  to  those  made  in  Section  III 
in  deriving  "optimal"  receiver  structures: 
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(1) 

(2) 

(3) 

(4) 

As  previously  noted,  these  assumptions  depart  somewhat  from  reality, 

but  will  serve  for  the  present  heuristic  analysis. 

N 

We  also  assume  that  the  data  sequence  {d^}  ^  consists  of  independent 
binary  values,  being  +1  or  -1  with  equal  probability;  and  that  the  channel 
noise  is  independent  of  the  multipath  var iib1 as. 

B.  The  Variance  of  7,(0) 

Attacking  the  signal  component  v(0)  in  (43)  first,  suppose  that  the 
desired  datum  in  v(0)  —  i.e.,  d^  —  is  +J .  Then,  using  the  assumptions 
discussed  above,  it  is  straightforward  to  show  that 


K—l  K~1  K— 1 

{a,  ,  {0.  }.  and  { t,  are  independent  sets  of  random 

k  0  k  0  k  0 

variables  and  all  variables  In  such  set  are  independent  of 
each  other. 

K—l  K—l 

The  number  of  paths  K  is  independent  of  (a^Jg  »  {Qj^g  and 

{t  J*"1 
1  k*0  ' 

All  t.  's  are  equi-distributcd  over  [0,A]  with  probability 

K 

density  p(tjc);  see  the  discussion  in  Section  IIIC. 

Each  0^  is  uniformly  distributed  over  (0 , 2tt ] . 


E[v(0) 


d  =+l 
o 


Y 


(0)ek[ 


K-l 

E 


k=0 


2£a2 


(44) 


Here,  an  overbar  denotes  "expectation,"  and  E^,  denotes  expectation  over 
the  random  variable  K.  We  have  used  the  fact  that  y(0)  =  2^,  being 

the  energy  in  the  transmitted  signal  s(t)  =  Re[o(t)e^Ut] ;  we  have  also  let 


,  .  K-l  - 

i“J2i 

k=0 


(45) 
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Not  so  straightforwardly ,  one  can  show  (see  Appendix  A)  that,  when  WA  >>  1, 


var (v(0) J 


-CO  ^ 

=  (2£)2B2+[|  |S(v)  |4dv)  [  53  q(nT)]C2 

=+l  n=-N 


where 


o  K-l  , 

,2  „  ,  i  i  4 


K-l 


b"  =  uE  Kl  -  E  l\l  i 

*  k=0  K  k=0 


K-i 

=  ek(£  var ( lakr ) i 


k=0 


9  K-l  K-l  - — - —  K-l  - —  2  K-l  - 

c  -  EKtE  E  l“kl  l«J  ]  =  ek[(E  l“kl  )  -  E  l\l 

k=0  £=0  K  k=0  k=0  k 

kjtJl 


and 


j: 


q(t)  =  1  p(t)  p(T-t)dt. 
0 


(46) 


(47) 

(48) 


(49) 


The  variance  in  (46)  is,  in  fact,  independent  of  the  condition 
dQ  =  +1 .  It  is  important  to  note  that  the  two  contributions  to  this 
variance  have  distinctly  different  interpretations.  The  first  term  is 
independent  of  N  and  would  be  present  even  if  there  were  no  multipath- 
induced  interference;  it  reflects  the  variation  of  v(0)  around  its  mean 
due  to  the  fluctuations  in  the  strengths  of  the  paths  that  have  been 
combined  by  the  RAKE  receiver  to  determine  the  datum  dQ  in  v(0).  The 
second  term  is  due  to  multipath- induced  interference;  it  reflects  the 
fluctuations  in  v(0)  due  to  the  presence  of  both  interpath  interference 
on  the  "present"  transmission  (n=0)  and  intersymbol  interference  from 
neighboring  transmissions  (n^O) .  The  second  term  will  vanish  only  if, 
with  probability  one,  there  is  but  one  path  (K=l)  in  (48). 

We  next  consider  the  term  n(0)  in  (4J).  It  has  zero  mean,  so  its 
mean  adds  nothing  to  (44).  Its  variance  is  derived  in  Appendix  B.  Since 
v(0)  and  n(0)  are  assumed  to  be  independent,  the  variance  of  z(0)  becomes 


ifitl 
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var  [  z (0) ]  =  var[v(0)]  +  var[n(0)] 

Joo  N 

|S(v)j4dv][  £  q(nT)]C2  +  4£noA2,  (50) 

— ®  n=-N 

where  N  is  the  (single-ended)  power  density  of  the  channel  noise,  assumed 
o 

to  be  white. 


C .  A  Measure  of  Performance  Degradation 

As  mentioned,  only  the  second  term  in  (50)  is  due  to  the  multipath- 
induced  interference  of  interest  to  us  here.  In  the  next  section,  for 
the  purpose  of  estimating  error  probabilities,  we  shall  model  this 
interference  as  approximately  Gaussian,  equivalent  to  an  additional 
amount  of  channel  noise.  In  this  framework,  multipath-induced 
interference  degrades  system  performance  by  the  same  amount  as  would  an 
increase  in  channel  noise  power  by  a  factor  of 


D 


1  + 


|S(v)|4dv][  £  q (nT) ) 
n=-N 


(51) 


We  now  specialize  (51)  with  some  additional  assumptions.  First, 
suppose  that  |s(f)|  is  flat  over  the  band  |f|  ^  and  zero  elsewhere. 

Then,  since  the  integral  of  |s(f)|2  is  equal  to  2^t, 


!  s(f )  | 2 


and 


f 


at> 


—co 


|s(f) 


(52) 


(53) 


Second,  suppose  that  A  >>  T,  so  that  we  can  use  the  following  approximation 
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where  £  =  |a|^£is  the  mean  energy  per  bit  arriving  per  path.  Thus, 

the  degradation  due  to  multipath-induced  interference  increases  with 

k£/Nq  (as  this  interference  increasing  dominates  the  channel  noise)  and 

decreases  with  TW  (as  the  sidelohes  of  y(t),  hence  the  interference 

contribution  by  each  interfering  path,  decrease). 

As  an  example,  suppose  —  using  parameters  typical  of  the  simulation 

£ 

experiments  in  Section  VI  —  that  K  5  25,  TW  =  127  and  —  =  5  dB.  Then 

N 

o 

D  *  1.62,  or  2.1  dB,  not  much  degradation  to  pay  for  the  ability  to 
increase  the  data  rate  by,  say,  a  factor  of  10! 

We  close  this  section  by  noting  that  if  the.  errors  in  the  estimates 
{a^},  i ! ,  lt^}  are  taken  into  account,  then  —  as  shown  in  Appendix  C  — 
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another  term  having  value  T  is  added  to  (57),  where  ^  is  the 

sounding  signal's  energy.  By  ignoring  estimation  noise,  we  have  implicitly 

assumed  that  £  /S  >>  A/T. 
x 

Although  tne  analysis  above  was  performed  for  a  purely  coherent  PSK 
system,  we  shall  assume  in  Section  VI  that  the  main  results,  (51)  and 
(57),  carry  over  approximately  to  a  DPSK  system. 


Exact  evaluation  of  the  error  probability  of  a  multipath  receiver 
quickly  becomes  intractable  as  the  assumed  model  for  the  multipath 
channel  becomes  more  realistic.  Most  analyses  assume  independent, 

Rayleigh-d is tribute  ^ath. strengths,  although  some  also  involve 
Rice  [13]  or  Nakagami  [  4  ]  strength  distributions  and/or  correlated 
strengths  [21 1.  [31  ] .  With  some  exceptions  [  4  ] ,  [31],  the  path  strengths 
are  assumed  to  be  equi-distr ibuted  and  the  path  delays  known.  None  of  the 
cited  analyses  takes  into  account  the  multipath-induced  interference 
discussed  in  the  previous  section,  although  Monsen  [15],  [lb],  [17]  and 
Morgan  [18]  consider  such  interference  in  a  different  context.  Exact 
analytical  evaluation  of  the  error  probability  for  the  model  in  Section  II 
—  involving  unknown,  non-Poisson  path  delays;  correlated  Nakagami-distributed 
path  strengths;  spatial  correlations  and  i nhomogeno i ( les ;  and  multipath 
induced  interference  —  seems  beyond  tractability. 

Another  approach  is  to  evaluate  the  error  probability  of  various 
receivers  by  simulation  techniques,  using  the  channel  simulation  program 
described  in  Section  II;  the  results  of  some  such  simulation  experiments 
are  given  in  Section  VI.  In  the  present  section,  however,  we  shall 
present  a  much-simplified  analysis.  This  analysis  will  provide  both 
intuitive  insight  into  the  performance  of  several  receivers,  and  a  basis 
for  determining  (in  Section  VI)  the  gap  between  reality  and  the  results 
of  rather  standard  analytical  methods.  Tn  our  analysis,  we  assume  lack 
of  knowledge  of  path  delays,  but  still  use  very  much  over-simplified 
statistics  for  the  path  delays  and  strengths. 

We  limit  ourselves  here  to  a  system  using  DPSK  transmission  and 


simplifications  of  the  receivers  of  Figures  11  and  17.  The  analysis 
depends  heavily  upon  the  work  of  Charash  [3],  [4]. 


A  simplification  of  the  receiver  of  Fig.  11  is  shown  in  Fig.  21, 
where  the  pertinent  portion  of  Fig.  10  is  also  included.  Notice  that 
the  (generally  time-varying)  nonlinearity  F [ • , • ]  has  been  eliminated, 
i.e.,  replaced  by  a  time-invariant  linear  device.  For  the  case  of  Rayleigh 
path  strengths  with  equal  mean-square  values,  the  nature  of  the  approxi¬ 
mation  is  shown  in  Fig.  12  (see  also  (19)  and  (24)).  By  eliminating  F, 
we  are  giving  up  the  strong-path-accentuation/ noise-suppression  feature 
of  the  receiver  that  is  shown  by  the  waveforms  of  Fig.  13(c),  (d). 

In  Fig.  21,  we  also  have  taken  the  integrating  filter's  impulse 
response  h  (t)  to  be  of  the  form  of  (26).  If  this  h^(t)  were  optimal, 
it  would  mean  that  paths  occur  with  uniform  density  over  the  interval 
[0,A];  see  Section  IIIC,  where  we  would  have  p(t^)  =  1/A,  0  £  t^  £  A.  In 
actuality,  paths  do  not  occ  uniformly,  so  this  h^(t)  is  sub-optimal. 

Since  the  lowpass  filter  of  Fig.  10  serves  only  to  eliminate  the  double- 
carrier-f requency  terms  in  the  multiplier  output,  its  function  is 
handled  in  Fig.  21  by  hp(t),  also  a  lowpass  impulse  response. 

'Hie  decision  circuit  input  in  Fig.  21  will  still  look  much  like  the 
waveform  of  Fig.  13(e),  but  it  will  now  be  an  integration  of  y(t)  of 
Fig.  13(c),  rather  than  of  Fig.  13(d).  It  must  still  be  sampled  by  the 
decision  circuit  at  the  instants  shown  in  Fig.  13(f). 

The  simplified  sub-optimal  receiver  just  discussed  has  also  been 
considered  by  Fralick  [  7],  who  descriptively  calls  it  a  post-detection 
integrator  (PDI),  a  name  we  shall  use. 

The  error-probability  performance  of  the  PDI  can  be  estimated  from 
the  work  of  Charash  [31,  [4],  He  has  given  error-probability  expressions 
and  curves  for  a  system  having  the  following  characteristics: 


(1)  Binary  transmission  using  uniformly  orthogonal,  zero-sidelobe 


signals,  i.e.,  ones  for  which  |y^,,(t)j  of  (7)  is  identically  zero 
rather  than  merely  satisfying  (8a),  and  |y^(t)i,  i  =  1,2,  is 
identically  zero  for  |t|  >  1/W  rather  than  merely  satisfying  (8b). 

(2)  A  receiver  of  the  form  of  Fig.  7  (with  M=2),  but  with  the  non- 

-  2 

linearity  being  a  simple  time-invariant  square-law  device,  F(x,t)=x  . 

(3)  In  integrating-f liter  impulse  response,  h^(t),  of  the  form  of 
(26). 


Note  that  this  system  bears  the  same  relationship  Lo  our  DPSK/PDI.  system 
in  the  multipath  case  as  the  usual  orthogonally  keyed  binary  system  bears  to 
a  standard  DPSK.  system  in  the  single-path  case  (cf.  Figures  5  and  10). 

We  shall  exploit  this  relationship  to  estimate  the  error  probability  for 
the  DPSK/PDI  system. 

Charash  also  makes  the  following  assumptions  about  the  channel: 

(4)  All  paths  are  independently  and  identically  Nakagami- 
distributed. 

(5)  Exactly  K  paths  arrive  at  random  over  [0,A].  Of  course,  to 
preserve  resolvability  condition  (4),  we  must  have  K  <<  WA,  where 
W  is  the  transmission  bandwidth. 

(6)  A  <<  T,  so  there  is  no  ippreeiable  intersymbol  interference. 

(7)  The  additive  channel  noise  is  white  and  Gaussian,  with  single- 
ended  power  density  N  . 

Under  these  seven  assumptions,  Charash  has  derived  an  expression 

for  the  bit  error  probability,  Pp(K),  conditioned  on  the  number  of  paths 

K.  (The  unconditional  bit  error  probability  would  then  be  the  average 

of  P_(K)  over  the  Poisson  distribution  of  K  that  is  implied  by  (5).) 

E 

As  mentioned,  our  DPSK/PDI  system  is  related  to  Charash* s  system  in 
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the  same  way  as  an  ordinary  single-path  DPSK  system  is  related  to  an 
ordinary  single-path  orthogonally  keyed  system.  The  performances  of 
these  latter,  single-path  systems  are  well  known  [24]  to  be  exactly 
3  dB  separated,  with  the  DPSK  system  having  the  advantage;  i.e.,  the 
DPSK  system  requires  3  dB  less  SNR  to  achieve  the  same  bit  error 
probability  as  the  orthogonally  keyed  system.  We  now  conjecture  that  the 
same  3-dB  separation  holds  (at  least  approximately)  for  the  corresponding 
multipath  systems. 

Using  this  conjecture,  we  have  used  Charash's  results  to  estimate  the 
error  probability  of  the  DPSK/PDI,  of  course  still  assuming  that  assumptions 
(3)- (7)  hold  and  thrt  the  DPSK  signal  has  zero  sidelobes.  Figure  22 
shows  some  typical  results,  for  Rayleigh  fading  (a  special  case  of  Nakagami 
fading)  ana  for  WA  =  40.  The  curves  there  show  the  signal -to-noise  ratio 
£/Nq  required  per  path  (see  (57)  et  scq . )  in  order  to  achieve  a  given 
value  of  P  (K) ,  as  a  function  of  the  number  of  paths  K.  Note  that,  since 
typical  urban  multipath  spreads  are  A  =  4-5  ps,  WA  =  40  implies  a  system 
bandwidth  of  8-10  MHz;  and  since  we  postulate  that  A  <  -  T  for  the  PU1, 
we  are  assume  a  spectrum  spreading  factor  of  TW  >>  40. 

On  the  K=1  axis,  we  have  shown  for  reference  the  performance  of  the 

* 

optimal  single-path  DPSK  receiver  of  Fig.  10.  One  can  see  that  if  a 

multipath  medium  has  a  number  of  Rayleigh  paths  of  more-or-less  equal 

mean-square  strengths,  our  DPSK/PDI  system  will  do  substantially 

(5-20  dB)  better  than  a  receiver  that  locks  onto  a  single  path 

(say  the  LOS  path)  and  disregards  the  remaining  paths.  This  comparison 

strikingly  illustrates  the  diversity  performance  improvement  that 
_  __________ 


This  is  derived  by  averaging  the  nonfading,  single-path  DPSK  error  prob¬ 
ability  [24]  of  (1/2)  exp(-£/N^]  over  a  Rayleigh  distribution  for  /& 
having  rins  value  /'j,  to  obtain  P  =  [2(1  +  £VNq)]  ^  where  Q,  = 


multipath  channels  can  afford,  as  discussed  at  the  beginning  of 
Section  III. 

It  must  be  stressed  that  the  curves  of  Fig.  22  take  into  account 
neither  intersymbol  interference  nor  interpath  interference  (see  footnote, 
page  33).  The  former  is  precluded  by  the  assumption  A  <<  T,  the  latter 
by  the  assumption  of  a  zero-sidelobe  signal. 

B.  A  Digital  RAKE  Receiver  (DRAKE) 

The  channel-sounding  receivers  of  Figs.  16  and  17  involve  adjustment 
of  the  tap  gains  of  the  transversal  filter(s)  with  estimates  of  path 
strengths.  An  arrangement  that  is  much  simpler  to  implement,  although 
it  will  not  perform  as  well,  is  as  follows:  connect  each  tap  to  the 
summing  bus  or  not  according  as  it  is  estimated  that  a  path  is  present 
at  that  delay  or  not.  Only  a  transfer  of  the  estimates  (t^.)  of  the 
path  delays  is  then  needed  from  the  sounding  receiver  to  the  trans¬ 
versal  f liter (s).  We  call  this  simplified  configuration  a  digital  RAKE 
or  DRAKE  receiver. 

Figure  23  shows  a  DPSK/DRAKE  receiver  based  on  Fig.  17,  in 
which,  as  a  further  simplification,  the  square  rooter  has  been  eliminated. 
The  sounding  receiver's  estimates  of  the  path  delays  turn  the  associated 
tap  switches  on  the  transversal  filter  on,  connecting  the  taps  to  the 
summing  bus.  Note  that  if  all  tap  switches  were  turned  an,  Fig.  23  would 
revert  approximately  to  the  PDI  of  Fig.  21,  with  the  A-second  integration 
replaced  by  an  approximating  A-second  discrete-time  summation. 

In  analyzing  the  DPSK/DRAKE  receiver,  we  assume  that  the  sounding 
receiver's  estimates  t^  of  the  t^  are  exact  (t^t,),  and  that  these 
estimates  correspond  exactly  to  available  discrete  tap  delays  in  the 


transversal  filter.  The  receiver  then  takes  the  form  of  an  equal-weight 
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path  combiner  for  known  path  delays.  Charash  [3],  f A]  has  given  expressions 
and  curves  for  such  a  receiver  under  assumptions  (1),  (A),  (5),  (6)  and 
(7)  of  the  previous  section.  We  can  convert  these  results  to  the  DPSK 
format  of  interest  to  us  by  the  same  conjecture  as  used  previously,  i.e., 
that  Charash's  orthogonally  keyed  system  performs  3  dB  worse  than  our 
DPSK  system.  On  this  basis,  curves  of  the  required  £i/NQ  per  path 
necessary  to  achieve  a  given  P  are  shown  in  Pig.  2A  as  functions  of  the 
number  of  paths  K.  Only  curves  for  the  Rayleigh-fading  case  are  shown 
here;  these  could  also  have  been  obtained  from  [31]. 

As  in  the  case  of  Fig.  22,  the  curves  of  Fig.  2A  take  into  account 
neither  intersymbol  interference  nor  the  interpath  interference  in  the 
current  symbol.  Both  of  these  introduce  an  effective  additional  noise 
term.  The  performance  degradation  caused  by  this  term  has  been  estimated 
in  Section  IV  for  an  ideal,  coherent  PSK/RAKE  system.  In  the  next  section, 
we  shall  assume  that  his  same  degradation  holds  approximately  for 
DPSK/DRAKE. 
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VI .  Simulation  Experiments 

In  order  to  evaluate  the  error  probabilities  of  various  receivers  in 
the  "actual"  urban  multipath  environment  and  Lo  assess  the  applicability 
of  the  simplified  analyses  of  the  previous  two  sections,  a  number  of 
simulation  experiments  were  performed.  These  experiments  combined  the 
urban  propagation  program  SURP  described  in  Section  IID  and  various 
binary  modulator /demodulator  programs  described  below.  Two  sets  of 
experiments  were  performed:  "low  rate"  at  78.7  kbps  and  high  rate  at 
787  kbps.  Both  sets  were  restricted  to  DPSK  transmission  and  the  dense- 

ic 

high-rise  environment  (Area  A),  although  other  experiments  are  planned. 

A .  Multipath  Data  Base 

As  a  first  step,  a  data  base  of  urban/suburban  multipath  profiles  were 

kk 

obtained  by  using  SURP.  For  each  of  the  four  areas  A,  B,  C,  D defined  in 
Section  IIA,  150,000  profiles  were  generated  along  a  track  3000-ft  long 
(see  Fig.  3),  with  uniform  spacings  of  0.02  ft,  using  1,-band  (1280  MHz) 
program  parameters.  In  order  to  reduce  the  cost  of  the  runs,  new 
{ x ^ }  and  {a^}  strings  were  generated  only  at  every  fifth  point  —  i.e., 
at  0.1-ft  (0.12  wavelength)  spacings  —  and  used  again  for  the  subsequent 
four  points.  (This  approximation  is  quile  Justified,  since  spatial 
correlation  distances  for  {t^}  and  (a^)  are  tens  of  wavelengths.)  However, 
since  the  correlations  distance  for  the  {0^}  string  is  a  fraction  of  a 
wavelength,  a  new  {0^}  string  was  generated  at  every  point. 

Some  typical  {x^}  strings  and  a^  and  0^  values  are  shown  in  Fig.  25. 
One  sees  the  {x^}  strings  beginning  to  decorrelate  after  about  one 
★ 

The  communication  system  simulation  programs  were  written  and  run  by 
M.  Kamil.  They  are  part  of  a  Ph.D.  thesis  in  progress,  the  aim  of  which 
is  to  explore  urban  mobile  digital  systems  by  simulation,  in  much  more 
depth  than  the  results  given  here. 

kk 

This  data  base  was  established  by  H.  Hashcmi . 
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wavelength  and  showing  substantial  decorrelation  after  10  wavelengths. 

Bin  13  happens  to  have  a  path  in  all  {x^}  strings  shown,  and  the 
corresponding  a^  and  6^  values  generated  for  that  bin  are  displayed. 

Again,  for  a^,  one  sees  partial  decorrelation  at  one  wavelength  and 
substantial  decorrelation  after  10  wavelengths.  On  the  other  hand  0^ 
begins  to  decorrelate  at  one-tenth  wavelength  and  shows  substantial 
decorrelation  at  one  wavelength. 

The  generated  profiles  were  stored  on  tape.  In  addition,  for 
future  use  in  planned  narrowband  experiments,  the  phasor  sums 

j  q  (tt)  j  ^ 

£a^n^e  k  _  n  were  computed  and  stored  for  each  n  (n  =  1 ,  .  .  . ,  150, 000) . 
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B.  "Low-Rate"  Experiments 

Recall  from  Fig.  2  that  the  simulated  multipath  profiles  are  discrete¬ 
time  impulse  responses  of  a  lowpass-equivalent  channel  filter,  with  time 
running  in  100-ns  increments.  In  designing  the  "low-rate"  communication 
system  simulation  experiments,  we  postulated  a  DPSK  transmitted  signal 

whose  basic  lowpass  equivalent  waveform  s(t)  is  a  127-chip  maximal -length 

* 


shift  register  (MLSR)  sequence  having  a  100-ns  chip  duration. 

126 

the  chip  sequence  by  {s.J^  >  the  lowpass-equivalent  transmitted 

for  an  N-bit  transmission  is  then  of  the  form 


Denoting 

signal 


“i  ■  £„  Vl-127„-  1  -  0 . 127N'1- 

n=0 


(58) 


Here,  s,  is  defined  to  be  zero  for  j  outside  [0,126],  and  d  d  ,  =  +1 
j  n  n-1  — 

according  as  bit  n  (n=l,...,N)  is  0  or  1,  respectively  (see  Fig.  13(a),  (b)) 


The  use  of  an  MLSR  sequence  has  a  beneficial  effect  in  our  context: 
whenever  s(t)  is  sent  twice  with  the  same  polarity,  the  right-hand  sidelobe 
fluctuations  of  the  matched  filter's  first  response  cancel  the  left- 
hand  sidelobe  fluctuations  of  the  second  response,  thus  considerably  reducing 
interpath  and  intersymbol  interference. 
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The  sequence  was  convolved  with  a  space-varying  sequence  of 

multipath  profiles  to  form  a  channel  output  sequence  that  simulated  the 

effect  of  vehicle  motion.  To  this  were  added  independent  complex-valued 

2 

Gaussian  noise  samples  having  variance  o^,  which  was  made  a  program 
parameter  to  allow  adjustment  of  the  SNR.  The  noisy  received  sequence  was 

1  o  A  A 

then  convolved  with  the  matched  filter  impulse  response  {s.„,  . } 

126-i  i=0 

At  this  stage,  one  has  a  discrete-time  representation  of  the 
matched-filter  output  in  Fig.  21  or  Fig.  23.  This  output  was  then  DPSK- 
demodulated  by  multiplying  it  by  the  complex-conjugate  of  a  127-chip- 
delayed  version  of  the  same  output  and  taking  the  real  part  of  the  product. 
No  nonlinear  enhancement  (e.g.,  as  in  Fig.  11)  was  attempted.  The  DPSK- 
demodulator  output,  which  is  a  noisy, discrete-time  version  of  Fig.  13(c), 
was  processed  directly  by  various  decision  mechanisms,  to  be  described. 

The  signal  format  wo  are  using  in  this  "low  rate"  case  implies  the 
following  parameters:  signal  bandwidth,  W  =  (100  ns)  1  =  10  MHz;  spectrum 
spread  factor,  TV!  =  127;  bit  rate,  1/T  =  78.7  kbps.  Further,  since  the 
multipath  profiles  are  restricted  to  71  samples  spaced  by  100  ns  (see 
Fig.  2),  we  have  WA  =  70.  Thus,  A/T  =  0.35,  so  that  intersymbol  inter¬ 
ference  is  negligible. 

Notice  that  a  vehicle  moving  even  at  60  mph  will  move  only  about 

_3 

0.001  ft  (about  10  wavelength)  per  bit,  so  bit-to-bit  changes  in  the 
phases  {0^}  negligibly  affect  the  DPSK  demodulation.  We  therefore 
simplified  the  experiments  by  "sampling"  the  transmissions  periodically 
in  space.  That  is,  we  imagined  that  the  vehicle  listened  to  10  bits  at 

•k 

In  fact,  for  computational  convenience,  the  order  of  the  operations  was 

quite  different:  the  convolution  of  { s , }  with  {s,„.  ,)was  done  first, 

1  1/b-i 

then  the  modulation  with  the  dn  sequence,  then  the  convolution  with  the 
multipath  profile,  then  the  addition  of  appropriate  noise. 
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each  of  10,000  spatial  positions  spaced  0.3  ft  (0.39  wavelength)  apart, 

for  a  total  of  100,000  bits  per  experiment.  In  this  way,  we  obtained 

-4  -5 

adequate  data  to  determine  bit  error  probabilities  down  to  10  -  10  , 

while  guaranteeing  a  representative  geographical  cross-section  (3000  feet) 
of  multipath  responses  and  a  representative  ensemble  of  noise  waveforms. 

Seven  decision  mechanisms,  operating  on  the  simulated  DPSK  demodulator 
output,  were  evaluated.  These  are  described  below,  and  should  be  visualized 
in  connection  with  the  waveform  y(t)  in  Fig.  13(c). 

(1)  First  path  (FP):  this  is  a  classical  mechanism,  in  which,  in 
each  bit  interval,  the  first  threshold  crossing  of  |y(t)|  is  examined, 
and  a  0  or  1  is  decided  according  as  y(t)  is  positive  or  negative. 

In  our  experiment,  the  threshold  level  was  set  at  a  standard  deviation 
of  the  matched-filter  output  noise,  and  we  idealized  the  mechanism  by 
assuming  that  the  first  crossing  was  always  due  to  a  signal  peak. 

The  position  of  the  first  crossing  was  found  by  having  the  program 
successively  compare  a^.a^ja^,...  with  a  subsidiary  threshold  that 
was  proportional  to  the  threshold  on  |y(t)|.  Because  of  this 
idealization,  our  simulation  result  is  a  lower  bound  on  the  per¬ 
formance  of  the  actual  mechanism. 

(2)  Largest  path  (LP) :  Here,  the  largest  peak  in  |y(t)|  in  each 
bit  interval  is  examined  and  the  polarity  of  this  peak  determines 
the  output.  Again  we  idealized  (not  much  of  an  idealization  in  this 
case)  by  assuming  that  the  largest  peak  is  a  signal  peak.  Its 
position  was  identified  as  that  of  the  largest  a^. 

(3)  Post-detection  integrator  (PDI) :  This  is  the  decision  mechanism 
of  Fig.  21.  However,  instead  of  integrating  over  a  full  A  =  7  us 
interval,  we  integrated  only  over  a  5-us  interval  beyond  the  LOS 


path. 


78 


(4)  Adaptive  post-detection  integrator  (APDI) :  Here,  the  PDI 
integrator  was  made  adaptive  by  letting  the  integration  interval 
run  only  from  the  first  crossing  to  the  last  crossing  of  |y(t)| 
above  a  threshold.  The  threshold  level  and  the  method  of  determining 
the  first-  and  last-crossing  positions  were  identical  to  those  used 
for  the  FP  mechanism. 

(5)  Weighted  post-detector  integrator  (WPDI):  Recognizing  that 
later  paths  in  a  multipath  profile  are  likely  to  be  weaker  and  less 
frequent  than  earlier  paths,  we  here  changed  the  impulse  response  of 
the  integrator  h^O)  of  Fig.  21  from  a  box  car  (equation  (26))  to  a 
decreasing  exponential  with  a  1.5-ps  time-constant  (an  RC  integrator). 

(6)  DRAKE:  This  is  the  receiver  of  Fig.  23.  The  switch  positions 
were  determined  by  assuming  that  the  sounding  receiver  will  identify 
as  paths  only  those  that  are  strong  enough  to  cause  |y(t)|  in  the 
data  receiver  to  exceed  a  threshold.  The  threshold  level  and  the 
method  for  identifying  above-threshold  paths  are  as  in  the  FP 
mechanism. 

(7)  RAKE:  Here,  Fig.  23  is  modified  to  restore  the  ideal  tap 
amplifiers  of  Fig.  14  in  place  of  the  switches.  The  amplifier 
gains  were  determined  by  assuming  that  the  sounding  receiver's 
estimates  of  the  a^'s  are  perfect.  (See  Appendix  C.) 

Mechanisms  (1),  (4)  and  (6)  depend  on  there  being  at  least  one  path 

above  threshold.  At  very  small  signal-to-noise  ratios,  there  were  a 
small  number  of  multipath  profiles  for  which  no  path  exceeded  the  threshold, 
and  the  demodulators  deleted  the  corresponding  sequence  of  ten  bits. 

Since  these  bits  by  definition  are  very  noisy,  ignoring  them  leads  to  an 
optimistic  estimate  of  the  error  probability,  P^ .  The  effect  was  not 
serious,  however,  causing  less  than  a  10%  bias  for  P  >  0.1  and  no 
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noticeable  bias  for  <  0.1.  But,  in  order  to  assure  valid  comparison 
of  mechanisms  (1),  (A)  and  (6)  with  (2),  (3),  (5)  and  (7),  which  have 
no  threshold  but  which  would  be  adversely  affected  by  inclusion  of  the 
noisy  bits,  we  deleted  these  bits  from  the  latter  demodulator  outputs 
also . 

The  empirical  bit  error  probability  curves  determined  from  the 

simulation  are  shown  in  Fie.  26  as  functions  of  both  fc  /N  and  £,„„/N  • 

av  o  LOS  o 

Here,  Nq  is  the  (single-ended)  channel  noise  power  density,  as  defined 

in  Section  IV.  £  is  the  average  energy  received  per  bit  per  path, 

determined  by  calculating  the  mean-square  path  strength  of  all  paths  in 

all  profiles  used  in  the  experiment.  t*ie  avera8e  energy  received 

per  bit  on  only  LOS  paths,  determined  by  calculating  the  mean-square  LOS 

path  strength  in  all  profiles  having  an  LOS  path.  In  actual  system 

design,  is  a  quantity  derivable  from  the  power  budget,  being 

determined  from  free-space  and  excess-attenuation  calculations  [33]. 

For  comparison,  also  shown  in  Fig.  26  are  classical  error-probability 

curves  for  an  optimal  DPSK  system  operating  through  a  single-path  channel, 

either  nonfading  or  Rayleigh  fading  (see  footnote,  page  67).  For  these 

curves,  the  abscissa  is  where  is  the  average  energy 

received  per  bit  over  the  single  path:  the  /N  scale  has  no  meaning 

i  n  '“av  o 

for  the  comparison  curves. 

From  Fig.  26,  we  can  draw  the  following  conclusions: 

(1)  Comparison  of  the  theoretical  one-path  Rayleigh-fading  curve 
with  the  empirical  f irst-path-above-thrcshold  (FP)  curve  shows  that 
patli  fading  (at  least  of  initial  paths)  Is  worse  than  Rayleigh,  a 
conclusion  already  drawn  by  Suzuki  [26]  and  Charash  [4].  (The 
fact  that  the  FP  curve  lies  below  the  Rayleigh  curve  for  small  SNR 


merely  shows  the  value  of  searching  for  an  above-threshold  path  when 


SIGNAL- TO-NOISF.  RATIO  (dB) 


Figure  26.  Krror  probabil  ities  for  the  "low  rate"  ease,  from  simulations 
(TW-1 2 7 ,  A/T-0.55,  R=l/T=78.7  khns.J 
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the  noise  is  substantial,  rather  than  always  utilizing  the  same 
path.)  Comparing  of  the  FP  curve  witli  Charash's  single-path  P 

E 

curves  [4]  —  as  appropriately  modified  for  DPSK  signalling  — 
shows  that  the  F1’  curve  reflects  Nakagami  fading  with  parameter 
m  =  0.75,  rather  than  m=l  (Rayleigh  fading). 

(2)  Comparison  of  the  WPOI  and  PDI  shows  that  a  h_ps  box-car 
integrator  out-performs  a  1.5-ps  RC  integrator,  at  least  for  dense 
urban  multipath.  However,  there  is  reason  to  believe  that 

other  urban  areas  and  for  other  RC  time  constants,  this  conclusion 
might  be  reversed. 

(3)  An  LP  receiver,  which  "captures"  the  largest  peak  in  y(t), 

works  overwhelmingly  better  than  one  that  captures  the  first 

available  peak.  (This  is  the  familiar  "largest  of"  voting  diversity 

effect.)  In  fact,  the  FI’  receiver  works  less  than  one  dB  worse 

-5  -2 

than  a  PDi  in  the  range  10  _<  P  £  10  ,  provided  that  the  largest 

peak  in  y(t)  indeed  corresponds  to  the  largest  path  rather  than  to 
noise.  However,  we  note  that  the  PDI  is  easier  to  implement  than 
the  FP  receiver,  among  other  things  requiring  less  accurate  sampling 
and  timing  circuitry. 

(4)  Progressing  through  successively  more  complex  systems  (PDI, 

APDI,  DRAKE  and  RAKE)  gains  oi  the  order  of  one  dll  per  stop,  provided 
that  Al’Dl  ,  DRAKE  and  RAKE  have  perfect  estimates  of  the  path 
information  they  need. 

(5)  Because  of  their  path-combining  characteristics,  by  which  they 
accumulate  the  energies  of  many  almost-indepcndent  paths,  the  WPD1, 
PDI,  Al’Dl,  DRAKE  and  RAKE  receivers  operate  at  an  enormous  advantage 
over  a  receiver  that  utilizes  only  a  fixed  single  fading  path.  At 
small  SDR,  they  even  moderately  outperform  a  receiver  that  utilizes  a 
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single  nonfading  path;  while  they  lose  this  moderate  advantage  at 
large  SNR,  their  performance  still  comes  within  a  few  dB  of  the 
nonfading  single-path  case.  (See  Monsen  [15]  for  a  similar  result 
obtained  by  equalization  techniques.) 

If  equipment  simplicity  is  taken  into  account  as  well  as  performance, 
it  would  seem  desirable  to  choose  the  PDI  receiver  in  this  "low  rate" 
case,  for  at  the  expense  of  less  than  2  or  .3  dB  in  SNR  the  PDI  can  be 
substituted  for  the  vastly  more  complicated  DRAKE  or  RAKE.  APDI,  while 
relatively  simple,  outperforms  PDI  by  less  than  1  dB,  scarcely  worth  the 
additional  threshold,  AGC  and  integrate-and-dump  features  of  APDI. 

As  we  have  seen,  however,  the  PDI  structure  is  quite  susceptible  to 
intersymbol  interference.  Thus,  in  the  high-rate  case  (A>>T),  we  must 
rely  on  RAKE  and  DRAKE, with  all  their  complexities.  We  investigate  their 
high-rate  performance  in  the  following  section. 


C .  High-Rate  Experiments 

In  the  high-rate  experiments,  we  increased  the  data  rate  by  a  factor 
of  ten,  i.e.,  to  1/T  =  787  kbps.  Since  A  remains  at  7  ps,  we  now  have 
A/T  =  5.5,  so  considerable  intersymbol  interference  exists:  six  pre¬ 
decessor  and  six  successor  symbols  affect  the  present  symbol. 

We  wish  to  maintain  the  same  spectrum  spreading  factor,  TW  =  127, 
which  means  that  W  is  now  100  MHz.  (This  large  bandwidth  does  not 
require  exclusive  assignment,  however,  as  many  systems  with  different 
codes  could  co-exist  in  the  same  band.)  A  problem  is  posed  by  this 

increased  bandwidth,  for  it  implies  a  resolution  of  1/W  =  10  ns  in  the 

* 

Much  of  this  small-SNR  advantage  over  nonfading,  single-path  operation 
is  real,  reflecting  the  benefits  of  path  combination  when  the  noise  is 
large.  Some  of  the  advantage,  however,  reflects  the  previously  mentioned 
bias  in  the  experiment  at  small  SNR  due  to  discarding  certain  "bad" 
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multipath  model,  while  our  fundamental  data  and  channel  simulation 

program  are  based  on  100-ns  resolution.  Unfortunately,  virtually  no 

10-ns-resolution  data  on  the  urban  channel  at  the  frequencies  under 

consideration  are  available.  (See  f 10] .  however.) 

We  have  resolved  this  dilemma  by  an  approximate  extrapolation  of  our 

100-ns  data.  Each  100-ns  bin  in  Fig.  2(a)  was  divided  into  ten  sub-bins. 

If  Tj  =  1  for  some  if 0,  the  corresponding  path  was  placed  at  random  in 

one  of  the  ten  sub-bins;  and  if  x  =  1  --  i.e.,  an  I.OS  path  exists  — 

o 

the  path  was  always  placed  in  the  first  sub-bin.  Thus,  implicitly,  we 
have  structured  the  multipath  to  preclude  an  average  of  more  than  one 
path  every  100  ns.  While  this  restriction  almost  certainly  varies  from 
reality,  it  allows  us  to  use  SURF  and  the  data  base  generated  by  it  for  the 
high-rate  case.  Hopefully,  the  resulting  simulation  results  will  not  be 

■k 

too  far  off. 

In  other  respects,  the  high-rate  experiment  followed  the  low-rate 
experiment.  DPSK  signaling  was  postulated,  using  the  same  MLSR  sequence, 
now  having  a  10-ns  chip  duration  and  a  1.27-ps  bit  duration.  The  multipath 
impulse  response  with  which  the  transmitted  signal  (58)  is  convolved  is 
now  710  samples  long  rather  than  71,  but  consists  of  more  than  90%  zeros. 
Only  RAKE  and  DRAKE  were  evaluated,  for  —  as  seen  in  Sections  IT  IF  and 
IV  —  only  they  are  capable  of  suppressing  intersymbol  interference. 

(WPDI,  PDI  arid  APDI  irretrievably  scramble  the  responses  to  different 
symbols,  while  FP  and  LP  present  substantial  difficulties  in  deciding 
which  is  the  first  or  largest  path  of  the  current  symbol  amidst  the 
interlaced  responses  to  several  symbols.) 


In  fact,  since  we  are  dealing  with  an  average  of  about  25  paths  per 
profile,  it  is  unlikely  that  use  of  a  refined  10-ns  model  —  having,  say, 
50-75  paths  per  profile  but  delivering  the  same  total  energy  —  would 
result  In  a  substantially  different  error  rate. 
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The  results  of  the  high-rate  experiment,  again  run  for  approximately 
100,000  bits,  are  shown  in  Fig.  27,  where  the  RAKE  and  DRAKE  low-rate 

results  are  also  shown  for  comparison. 

We  see  that,  as  predicted  in  Section  IV,  the  presence  of  substantial 

intersymbol  interference  does  not  drastically  curtail  system  operation. 

In  fact,  the  degradation  formula  (57),  although  derived  for  coherent 

RAKE  and  a  much  simplified  multipath  model,  proves  to  be  remarkably 

accurate.  The  average  number  of  paths  per  profile  encountered  in  the 

simulation  run  was  23.2.  If  we  allow  the  simplified  multipath  model  of 

Section  IV  to  have  K  =  23.2  paths  and  the  £/NQ  parameter  of  that  model 

to  be  identified  with  A  /N  of  Fig.  27,  then  —  setting  TW  =  127  —  D 

wav  o 

of  (57)  is  1.3  dB  at  £  /N  =3  dB,  2.0  dB  at  A  /N  =5  dB  and  2.8  d3 

av  o  ^av  o 

at  >-  /N  =  7  dB.  These  are  almost  exactly  the  gaps  between  the 
'-'av  o  J 

two  RAKE  curves  shown  in  Fig.  27. 

D .  Comparison  with  the  .S' imp  1  i  f  fed  Theory 

In  Section  V,  we  gave  theoretical  error-probability  curves  for  the 
low-rate  PDI  and  DRAKE  demodulators;  see  Figs.  22  and  24.  Since  these 
curves  were  based  on  a  number  of  very  simple  but  often-used  assumptions 
(a  constant  number  of  uniformly  distributed  Poisson  path  arrivals; 
independent  Rayleigh-distributed  path  strengths;  no  inte.rsymbol  inter¬ 
ference;  a  zcro-sidelobe  DPSK  signal,  i.e.,  no  Interpath  interference), 
it  is  of  value  to  compare  them  with  reality.  The  comparison  is  shown  in 
Fig.  28.  There,  we  have  redrawn  the  simulation-derived  low-rate  PDI  and 
DRAKE  curves  from  Fig.  26  and  the  high-rate  DRAKE  curve  from  Fig.  27. 

Also  shown  are  theoretical  low-rate  PDI  and  DRAKE  curves  taken  from  Figs. 

22  and  24  (as  extended  from  [  3]),  which  wire,  obtained  by 

extrapolating  Figs.  22  and  24  to  a  value  of  K  equal  to  the  average  number 

of  paths  per  profile  (23.7)  encountered  in  the  simulation.  The  theoretical 
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SIGNAL-TO-NOISE  RATIO  (dB) 


Error  probabilities  for  the  hi^h-rate  case,  from  simulations 
(TW=127 ,  A/T~5 . 5 ,  R=l/T=787  kbps.) 


high-rate  DRAKE  curve  was  obtained  by  assuming  that  (57)  --  derived  for 
ideal  coherent  RAKE  —  applies  approximately  also  to  DRAKE;  in  (57),  we 
used  K  =  23.7. 

We  see  that  the  simplified  theory  leads  to  results  that  are  between 

-4 

3  and  7  dB  optimistic  in  neighborhood  of  P  --  10  ,  and  less  so  at  larger 
error  probabilities.  This  optimism  is  due  to  the  combined  effects  of  the 
following  factors: 

(1)  As  we  have  noted  in  connection  with  the  FP  curves  of  Fig.  26, 
actual  urban  fading  is  worse  than  Rayleigh.  Reference  to  the  worse- 
than-Rayleigh  curves  of  [3]  show  that  the  Rayleigh  assumption  leads 
to  about  1-dB  worth  of  optimism. 

(2)  Actual  urban  paths  are  correlated,  rather  than  independent;  we 
estimate  that  assuming  independence  leads  to  another  1  d3  of  optimism. 

(3)  Because  of  spatial  inhomogeneities,  actual  path  strengths  have 
variances  that  are  smaller  in  local  areas  than  on  a  global  basis. 

The  simplified  theory  imagines  that  the  global  variances  apply  at 

each  local  site,  leading  to  a  larger-than-actual  "path  diversity"  gain. 

(4)  The  actual  DPSK  signal  used  did  not  have  zero  sidelobes,  so  more 

interpath  and  intersymbol  interference  exists  than  in  the  simplified 

theory.  We  estimate  from  Section  IV  that  ignoring  this  interference 

-4 

accounts  for  about  2.5  dB  of  optimism  at  =  10 

(5)  Finally,  since  the  theoretical  low-rate  DRAKE  curve  is  optimistic 
in  estimating  the  required  £av/NQ  for  a  given  Pg,  use  of  (57)  will 
lead  to  an  optimistic  value  of  D  to  be  applied  in  obtaining  the  high- 
rate  curve,  thus  compounding  the  error.  The  error  is  further 
compounded  by  the  assumption  the  (57)  applies  to  DRAKE,  since  Fig. 

27  shows  that  the  intersymbol  interference  is  more  damaging  to  DRAKE 


than  to  RAKE. 
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Thus,  the  simplified  theory,  although  based  on  common  analytical 

-2 

assumptions,  must  be  treated  with  caution,  especially  at  P_  <  10  ;  that 

is,  its  optimism  should  be  recognized  in  attempting  to  apply  it  to  urban 
multipath. 
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Appendix  A:  Derivation  of  (46) 

From  (42)  and  using  the  assumptions  listed  just  prior  to  that 
equation,  we  obtain 


N 


K-l 


E[  | v (0) ]  ] 


d  =+l  n=-N  k,2,q,r=0 
o 


=  E.  EK  .  X  „  VVVr  Y(nT-tk+tt)Y"(nT-tq+tr)  , 


(A-l) 


where  the  overbars  denote  expectation  and  K  denotes  expectation  over  the 
random  variable  K.  Note  that,  of  all  the  combinations  of  the  indices 
k,  2,  q  and  r  in  (A-l),  in  any  instance  where  one  is  distinct  from  the 


'h  ^ 

others,  the  expectation  a.  a. a  a  will  factor  into  the  form  a, a„a  a 

k  2  q  r  k  2  q  r 

(where  here  r  is  the  distinct  index);  this  factorization  follows  by  virtue 
of  the  independence  of  the  a^'s.  But,  because  of  the  assumed  uniform 
distribution  of  the  0  's  and  the  independence  of  the  a,  's  and  6  's, 

j0k 

*  Eta^e  1  =  0.  Thus,  the  only  surviving  terms  in  (A-l)  are  those  in 

which  indices  are  pairwise  equal,  and,  of  these,  those  where  k=r  and  2=q 
~  2  J20k 

vanish  because  =  Eta^e  J  =  Therefore, 


N  K-l 


K-l  K-l 


E[|v(0)|2]  =  £  Er[£  !«.  i  |y(nT)  \  +  £  £  la!"  htr  |Y(nT) 

n=-N  '  k=0  k=0  2=0 

k^S- 


+  E  E  l«kl2  K!2  M"T-tk+t£)|2, 

k=0  2= 0 
k/2 


(A- 2) 


where  we  have  noted  that  (A-l)  is  independent  of  dQ. 


Since 
Y(t)  »  ( 

J  —  o 


I  S(f )  I  2  eJ27Tftdf  , 


(A- 3) 


n>  *  >*  i  «  r,  -  — i 


we  have 
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l^nT-tj+tpl2  =  f  |S(f)j2|s(v)|2  oj27T(f_v)nT 

<  — oo 


j2n(f-v)tk  -j  2ti  (f-v) t^ 


(A-4) 


If  we  define  p(tk)  to  be  the  probability  density  distribution  of  the  t^'s, 


j2lTftk  fA  j2lTftk 

s  -  J  P(tk)  e  K  dtk  =  P(f) 


(A-5) 


is  the  characteristic  function  of  p(tk). 

2 

We  note  from  (23)  et  seq.  that  )y(nT)|  =  0  for  n/0.  Using  this  fact 

and  (A-4)  and  (A-5),  (A-2)  becomes 


E[  |  v  (0)  | 2  ]  =  (2g)2[E  £  |a  |4  +  F  £  £  |a  (2  |a  |2] 

k=0  k=0  1=0 


N  H  K-l  - j - o  Cf'“  9  ? 

E  E K  E  I  Kl  l«J  If  |s(f)|2|S(v)!2  |P(f-v) 
n=^N  K  k=0  H«=0  k  1  -® 


cj2i,(f-v)nT  dfdv> 


(A-6) 


where  we  have  used  the  fact  that  y(0)  =  2  a  £  being  the  energy  of  the 

transmitted  signal  s(t)  =  Re[o(t)e^wt] . 

From  (A-5)  we  see  that  for  a  reasonably  diffuse  p(t^)  —  in  the  sense 

of  Section  IIIC  —  the  "width"  of  P(f)  must  be  of  the  order  of  1/A. 

2 

On  the  other  hand,  the  width  of  |s(f)]  is  of  the  order  of  W.  But  we  are 

assuming  TW  >>  1  and  A  >  T,  so  W  >>  1/A.  The  situation  is  depicted  in 

2  2 

Fig.  29,  where  the  factors  |S(f)|  and  |P(f-v)|  in  the  integrand  of 


92 


(A- 6)  are  shown.  If  we  assume  that  |s(f)|  is  smoothly  varying, 
the  integral  on  f  in  (A-6)  van  be  written  approximately  as 


(  |s(f)  |2  |p(f-v)  |2  eJ27,fnTdf 

/  .O) 


=  is(v) | 2  f  |p(f-v) | 2 

j  CO 


J2r.fnT 
e  d  i 


|  S  (v)  (  2  q  (nT)  eJ2irvnT  t  jvj  _<  J 

I  i  W 

>  M  >  J 


00 

where  q(t)  =  (  |?(f) | 

J  — OU 


2  J  2irft., 
e  df 


j; 


p(t)  p(T-t)dT. 


The  double  integral  in  (A-6)  then  becomes 


(A-7) 


(A-8) 


q(nT) 


r  is<v)f 

J  — O 0 


dv 


(A-9) 


so  (A-6)  becomes,  after  a  little  manipulation. 
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'Ik1 
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(A-10) 


The  second  line  jn  (A-10)  is  recognizable  as  the  square  of  E[v(0)]  of  (44); 
subtracting  it,  we  get  var[v(0)],as  given  in  (46). 


■  f 
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Appendix  H:  Derivation  of  tin-  hast  Term  in  (50) 

j  27TfQt 

Suppose  that  the  channel  noise,  n(t)  =  Re[v(t)e  ] ,  is  a  bandpass 

white  process  having  spectral  density  N()/2  in  a  band  around  +f  .  The 
lowpass-equivalent  noise  v(t)  then  has  a  spectral  density  2N^  in  a  band  of 
the  same  width  around  1=0  (see  [31])-  v(t)  passes  through  the  receiver's 
cascade  of  the  matched  filter  and  transversal  filter,  which  has  overall 

■fe  /v 

transfer  function  S  (f)  H  ^(f)  —  see  (35).  Assuming  that  the  transversal 
filter's  taps  are  exactly  set  (a^  =  t^  =  t^,  all  k),  the  noise  output 

of  this  cascade  has  spectrum 


U(f)  = 


2N  |S(f) 


K-l  *  +J2*ftk 

L  \  e 

k=0 


9  K  K-l  *  j  2:i  f  ( t, -t  „  ) 

=  2N  |  S  (f )  j  2  £  £  Vl  e 

°  k=0  SL=0  * 


(B-l) 


Using  the  simplified  assumptions  on  multipath  statistics  listed  in 
Section  IV,  we  can  easily  show  that  the  average  of  U(f)  is 

U(f)  =  2No|S(f)|2  A2  ,  ( B— 2 ) 

where  A  is  as  in  (45).  The  variance  of  the  noise  at  the  output  of  the 
transversal  filter  is  therefore 

oo 

U(f)df  =  4£nq  a2  .  (B-3) 
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Appendix  C:  The  Effec t  of  Noi sv  Path  Ex t  imates 


In  the  analysis  of  Section  IV,  we  assumed  that  the  channel-sounding 

system  was  capable  of  supplying  perfect  estimates  of  path  variables  and 

that  the  RAKE  receiver  was  capable  of  using  these  exactly.  In  effect, 

we  therefore  assumed  that  the  sounding  receiver  can  supply  h  (t)  of  (30) 

or  H  (f)  of  (j4)  exactly  and  RAKE  can  then  configure  itself  into  a  filter 
m 

*  ~k  —  i  2  fi  f  A 

with  transfer  function  H  ^.(f)  =  H^(f)o  (see  (35)).  We  investigate 

the  implications  of  these  assumptions  here. 

Estimation  of  the  impulse,  response  of  an  unknown  medium  has  been 
investigated  in  [ 28  ] ,  where  the  optimum  soundi  .g  signal  and  optimum 
estimating  receiver  are  derived.  Generally  speaking,  the  optimal  signal's 
energy  is  distributed  in  a  way  depe.nding  on  the  noise  spectrum  and  average 
channel  transmission  function,  while  the  optimal  rc-eivur  uses  a  filter 
that  is  between  being  inverse  to  and  matched  to  the  optimum  signal. 

Mien  the  simplified  channel  model  of  Section  IV  is  assumed,  however, 
the  optimal  signal  distributes  its  energy  uniformly  over  the  transmission 
band  of  the  system  and  the  optimal  receiver  is  matched  to  this  signal. 

In  lowpass-equivalent  form,  the  sounding  receiver's  inatched-f liter  output 
is  a  noisy,  bandlimlted  version  of  (39),  i.e.. 


K-l 

k=0 


jC>k  sin  t, W ( t  —  t k ) 

V*  nw(T-T^)  +  i  (t)  ’ 


(0-1) 


where  t‘(t)  is  a  complex  lowpass  output  noise.  In  (C-l),  the  (sin  x)/>: 
pulses  have  been  scaled  to  have  the  same  areas  as  the  corresponding 
impulses  in  (30).  hence,  the  Fourier  transform  of  the  first  term  of 
(C-l)  agrees  with  that  of  (30)  within  the  band  |f|  W/2. 

We  approximate  the  RAKE  transversal  filter  by  assuming  that  its 
impulse  response  is 
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K-l 
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k“0 


*  sin  nW<A-t-tk) 
ak  TrW(A-t-"tk) 


+  5*(A-t) 


(C-2) 


instead  of  (33).  The  first  term  of  (C-2)  agrees  with  (33)  in  the  band 
|f |  <_  W/2  and  therefore  will  act  on  RAKE  input  signals  identically  as  (33); 
to  implement  it  would  require  infinitesimally  spaced  taps  on  the  RAKE 
delay  line,  however.  The  second  term  of  (C-2)  comprises  the  noise  in  the 
estimate  supplied  by  the  sounding  receiver. 

The  second  term  in  (C-2)  is  that  responsible  for  the  degradation 
in  performance  due  to  a  noisy  sounding  estimate.  As  mentioned,  the 
optimal  sounding  signal  for  our  case  has  a  flat  spectrum, 


|x(f)| 


|  f  |  1  W/2  , 


(C-3) 


where  is  the  sounding  energy.  In  order  to  have  the  scaling  of  (C-2), 
the  magnitude  of  the  frequency  characti ristic  of  the  sounding  receiver's 
matched  filter  must  be  /w/2fix  over  the  band.  The  noise  power  density 
of  £(f)  is  thus 


p{(f) 


fNo" 


1° 


|  f  I  1  w/2 

l  f |  >  W/2 


(C-4) 


since,  from  Appendix  B,  the  lowpass-equivalent  channel  noise  has  power 
density  2Nq. 

Let 


y(t)  “Cf3T"1(Hm(f)|S(f)|2ej2T,fT] 
K-l 

“  £  ak  Y(t-tk) 
k*0 


(C-5) 
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be  the  lowpass-equivalent  signal  output  of  the  data  receiver's  matched 

filter,  in  response  to  a  single  signal  transmission,  as  shifted  so  that 

its  LOS  peak  occurs  at  t=0.  The  RAKE  transversal  filter  convolves  this 

output  with  to  form  w(t),  the  signal  component  of  which  is 

specified  by  (36)  and  (38).  w(t)  also  has  a  component  due  to  noisy 

* 

estimation  of  channel  parameters,  obtained  by  convolving  the  £  (A-t)  of 

* 

(C-2)  with  y(t)  of  (C-5) .  This  component  of  w(t)  has  the  form 
t A  * 

1  £  (A-t)  y (t-r)dT  (C-6) 

•'O 

which,  when  shifted  to  the  left  by  A  seconds,  becomes 

t  &  * 

1  £  (x)  y(t+x)dr  .  (C-7) 

Jo 


Equation  (C-7)  gives  the  estimation-noise  component 
of  w(t)  when  only  a  single  signal  is  transmitted.  When  a  sequence  of 
signals  is  transmitted  with  PSK  modulation,  as  discussed  prior  to  (39), 
the  estimation-noise  component  becomes 


e(t)  =  1  £  (x)  £  d  y(t-nT+i)dr  , 

} 0  n=-N 


(C-8) 


where  we  have  indicated  that  the  responses  to  only  2N+1  symbols  are  in 
the  RAKE  delay  line  at  one  time. 

At  t=0,  when  the  signal  component  of  the  RAKE  filter  output  peaks 
for  the  present  symbol  (see  (42)),  the  estimation-noise  component  is 


e(0) 


N  .  A 

n^N  Mo 


£  (x)  y (x-nT)dx  . 


(C-9) 


There  is  also  a  term  in  w(t)  involving  convolution  of  £(t)  with  the  noise 
in  the  data  receiver's  matched-filter  output;  we  neglect  this. 
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This  has  zero  mean,  *.-<ce  (,  (t)  =  0,  and  variance 


N  N  -j -  - - - 

[|e(0)|2]  ■  V*  Z  dndm  \l  ^  (x)C(°)  y (x-nT)y  (o-mT)  dxdo 
n«-N  m--N  ™  0 


(C-10) 


Since  d^^  -  <$mn  (independent  data  symbols),  and  (from  (04)) 
- -  N  W 

C  (x)C(tf)  ■  -g-  S(t-o),  and  (from  (C-5)), 


K-l  K-l 


y (x-nT)y  (o-mT)  =  E  £  E  Vf  Y(x-nT-t.  )y  (a-mT-t.), 

k=0  i=0  £ 


(C-ll  ) 
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and  using  -  a^j  <5^  (from  the  simplified  multipath  model  of 


Section  IV),  we  have 


E[ |e(0) | 2]  = 


N  W  N  K-  l  /•  !- 

=  TT  E  EK  E  f 

Ox  n=-N  k=0-,l 


l«J  |y(i-nT-tk)j  dr. 


(C-12) 


Again  using  the  simplified  multipath  model  of  Section  IV,  in  which  all 
t^'s  are  independent  and  uniformly  distributed  over  [0,A],  we  have 


N  WA  N  ,A,A 


9  -  -  f  t  O 

E[|e(0)j  ]  =  ttt-  E  1  I  |Y(x-nT-t  )|2  drdt 
Ox  n*=-N  JQ'O 


(C-13 ) 


where  A  is  given  by  (45).  Using  the  transformation  u  =  T-t^,  v  =  T+tjc» 
whose  Jacobian  is  1/2,  (C-13)  becomes 


„  N  WA2  N  f  A  .  2A-  |u| 

E[  |  G  (0)  [  ]  “  E  \  dl  I  dv  |y  (u-nT)  J 

vOx  n*-N  * -A  *  ! u I 


(C-14 ) 
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«  nn  mi  i  - 

-p-7-  J]  l  (A- |u j )  | y (u-nT) |  du  . 
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Finally,  we  note  that  |y(t)|  is  a  narrow  pulse  of  width  1/W  <<  A. 
Treating  it  as  a  delta-function  of  area 
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|S(f)|4df 


(C-15) 


(since  the  Fourier  transform  of  y(t)  is  |s(f)|2,  the  signal's  lowpass- 


equivalent  energy  density  spectrum) , and  using  (53),  we  have 


„  4N  £2A2  N 

E[|e(0)l2]  ^  £  (A—  | n J T )  .  (C-16) 

Cx  n=-N 


When  T  >  A,  so  only  the  n**0  term  exists,  (C-16)  becomes 


2  4N  £2A2 

E[|e(0)!2]  =  — g -  .  (C-17) 


More  interestingly,  when  A  >>  T,  so  we  can  approximate  the  sum  in  (C-16) 

2 

with  an  integral  of  value  A  /T,  we  have 


E[ | e(0) [ 2] 


'-'X 


A 

T 


(C-18) 


2 

E [  | e (0) |  ]  is  another  term  to  be  added  to  (50),  adding  to  the 
degradation  caused  by  multipath-induced  interference  and  channel  noise. 
Taking  this  term  into  account,  the  increase  in  the  RAKE  output  noise 
variance,  as  given  by  (57),  now  becomes 


D  1  +  TW  N  P  T  ’ 
o  Ox 

where  we  have  assumed  that  A  >>  T  and  used  (C-18). 

We  thus  conclude  that  if  the  ratio  of  sounding-signal  energy  to 
data-signal  energy  satisfies 


(C-19  ) 
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(C-20) 


then  the  estimation-noise  contribution  to  the  total  RAKE  output  noise  is 


negligible.  If,  as  in  the  high-rate  case  of  Section  VIC 
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r  jr"  ’ 


A 

T 


then  £  should  be  10  times  or  more  greater  than  Q,  for  (C-20) 
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